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^ ABSTRACT 

o- 

' We present spectral energy distribution modelling of 6875 stars in uj Centauri, ob- 

taining stellar luminosities and temperatures by fitting literature photometry to state- 
of-the-art MARCS stellar models. By comparison to four different sets of isochrones, 
we provide a new distance estimate to the cluster of 4850 ± 200 (random error) ±1 20 
—i . (systematic error) pc, a reddening of E(B — V) — 0.08 ± 0.02 (random) ±0.02 (sys- 

■ tematic) mag and a differential reddening of AE(B — V) < 0.02 mag for an age of 12 
^sO [ Gyr. Several new post-early- AGB candidates are also found. Infra-red excesses of stars 

i were used to measure total mass-loss rates for individual stars down to ~ 7 x 10 -8 M Q 

■ yr" 1 . We find a total dust mass-loss rate from the cluster of 1.3 ±o - | xlO -9 Mq yr" 1 , 

with the total gas mass-loss rate being > 1.2 ±§| x 10" 6 M Q yr -1 . Half of the cluster's 

£NJ ■ dust production and 30% of its gas production comes from the two most extreme stars 

7— I ' - V6 and V42 - for which we present new Gemini/T-RcCS mid- infrared spectroscopy, 

possibly showing that V42 has carbon-rich dust. The cluster's dust temperatures are 
found to be typically ^ 550 K. Mass loss apparently does not vary significantly with 
metallicity within the cluster, but shows some correlation with barium enhancement, 
which appears to occur in cooler stars, and especially on the anomalous RGB. Lim- 
its to outflow velocities, dust-to-gas ratios for the dusty objects and the possibility 
of short-timcscale mass-loss variability are also discussed in the context of mass loss 
from low-metallicity stars. The ubiquity of dust around stars near the RGB-tip sug- 
gests significant dusty mass loss on the RGB; we estimate that typically 0.20-0.25 M 
of mass loss occurs on the RGB. From observational limits on intra-cluster material, 
we suggest the dust is being cleared on a timescale of ^ 10 5 years. 

Key words: globular clusters: individual: oj Centauri - stars: AGB and post-AGB - 
stars: mass loss - circumstellar matter - infrared: stars - stars: winds, outflows. 



1 INTRODUCTION 

1.1 Mass loss and globular clusters 

Galactic globular clusters (GCs) are unique stellar labora- 
tories, containing roughly co-eval populations of stars at 
kno wn distance s, covering the range [Fe/H] « -2.3 to so- 
lar (|Harrd H996V They are ideal locations in which to probe 
the later stages of development of low-mass (~0.8 Mq) stars, 
and their large (~ 10 5 — 10 6 ) populations mean they can 
harbour objects in very short phases of evolution. 

Stellar mass loss, and its implications for both stellar 
evolution and the cluster's fate are of particular impor- 



tance. All stars over ~0.8 Mq are thought to lose ^30% 
of their mass on the Red and Asymptotic Giant Branches 
(RGB/AGB) (e.g. 



Roodll 19731: iLee. Demarque fc Zinnlll994l : 



ICaloi fc D'Antona [2008!) . Warm giants may drive an out- 



flow of ~ 10 Mq yr via magnetic or hydrody- 
namic processes, (e.g. Dupree Hartmann. fc AvrettJ Il984l : 



M auas. Cacciari fc Pasauinil 12006). while the coolest, most 
luminous giants are thought to be able to sustain a wind 
of at least ~ 10 -6 Mq yr -1 through a combination of 
pulsa tions and radiation pressure on circumstella r dust 
(e.g. iGehrz fc Wooll Il97ll : iBowen fc Willsonl Il99ll h Mid- 
infrared (IR) spectra can be used to identify the species, 
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temperature, mass and nature of these circumstellar dust 
grains an d improve estimates of the total mass-loss ra te from 
the star ijLebzelter et al.ll2006l : | van Loon et al.ll2006h . 

Mass loss has important consequences for the remainder 
of the star's evolution and the chemical enrichment of inter- 
stellar material. Significant mass loss can decrease the num- 
ber of thermal pulses such a star undergoes, thus inhibiting 
the growth of the stellar core, and reducing dredge-up, which 
affects the chemistry and char acteristics of the post-AGB 
and planetary nebula phase (Vassi liadis fc Wood! 1 1993T ) . 
Mass-loss evolutionary processes also determine whether the 
star r emains an oxygen-rich AGB star or beco mes a carbon 
star (|van Loon et all Il998l ; iFrost et all Il998l ). Ultimately, 
mass loss will limit the white dwarf mass. Many aspects 
of the mass-loss process remain poorly understood. There is 
no clear consensus even on the primary driving mechanism 
behind mass loss, likely pulsation or continuum-driving of 
dust. By examining mass loss along the giant branches of 
populations of different metallicities, we can gain insight into 
the poorly-determined relationships between metallicity and 
the mass-loss rate, the wind speed and gas-to-dust ratio of 
the star s; and the correlation betwee n mass loss and grav- 
ity (e.g. iMcDonald fc van Loonll2007l ). The clearing of lost 
mass from within the cluster will also liberate mass from the 
cluster and could exacerbate its evaporation. Study of the 
clearing can also al low us to infer con ditions in the cluster's 
local environment (lOkada et al 1 l2007h . 

The Spitzer Space Telescope |Werner et alj 1 2004 
iGehrz et al.[|2007ft has, for the first time, allowed a com- 
plete census of dust production within GCs, due to its un- 
precedented sensitivity and angular resolution. A compre- 
hensive analysis of mass loss from stars in 47 Tuc was re- 
cently performed by lOriglia et al.l (|2007n (hereafter O+07). 
In this study, mid-IR flux excesses from Spitzer data were 
used to construct the mass-loss rate from individual stars 
in the cluster, attaining an empirical relation for mass loss 
from metal-poor giant stars. This relation suggests that sig- 
nificant dust production occurs along a considerable part of 
the giant branch. We compare our results to this analysis in 
Section 16.31 



1.2 Omega Centauri 

On the basis of our own Spitzer data, published in 
iBover et all lj2008h (hereafter B+08), and new mid-IR spec- 
tra of the two stars (V6 & V42) with the strongest IR excess, 
we here present an estimate of the mass-loss evolution in the 
most massive Galactic globular cluster: lo Centauri. 

u Cen is unique in the wealth of information avail- 
able - it is comparat ively nearby at a relatively well-known 



distance of ^5 kpc ([Ha rris 1996; van Lccuwcn et al. 2000; 
Ivan de Ven et al.ll2006l ; lDel Principe et alfeOOrJ); its high ra- 
dial velocity of vl s r = + 232 kms" 1 dMevlan fc Mavorll986l ; 
Ivan de Ven et ail 120061 ) allows easy spec troscopic mem- 
bersh ip and radial velocity determination (jvan Loon et al.l 
|2007| ; hereafter vL+07), and it also has independent mem- 
bership determinations fr om proper motion measurements 
l|van Leeuwen et al.l |2000| ; hereafter vL+00). The impor- 
tance of this is highlighted in B+08, which also contains an 
assessment of the effects of stellar blending on our longer- 
wavelength (lower-resolution) data. The cluster contains a 
statistically-significant sample of stars in most advanced 



stages of evolution (vL+07). Coupling this with a large exist- 
ing photometric dataset, it is possible to positively identify 
cluster members and probe the stellar population well down 
the giant branch at all wavelengths relevant to spectral en- 
ergy distribution (SED) modelling. 

The bulk of the stars in the cluster have a 

relatively low metallicity ■ [ Fe/H[ ~ -1.7 t o -1.6 

l|Norris. Freeman fc Mighelllll996l ; ISmith et alll200Ch . How- 
ever, a helium-e nriched, metal-intermediate popul a tion at 
[Fe/H l w -1.2 dNorris. Freeman fc Mighelll 1 19961 ; iNorrisI 
|2004| ) also exists; along with a metal-rich component, the 
'anomalous' RGB (RGB-a), with metallicities up to [Fe/H] 
~ - 0.7, which together comprise perhap s 10% of the clus- 
ter l|Lee et alj 119991 : iPancino et al.ll200ol . |2002h . Variations 
in surface abundances are also present, including oxygen- 
rich, M-type stars with titanium oxides; stars enhanced in 
CH and/or CN; and genuine carbon stars with molecular 
carbon (C2; vL+07). 

The origin of these sub-populations is undecided 
JSollima et alj 120051 : IStanford etall 120061 ; IVillanova et al.l 
120071 ) . though there is evidence that u> Cen may be the rem- 
nant nu cleus of a tidally disrupted dwarf sph eroidal galaxy 
(dSph) jZinnecker et al.lll988l ; lFreemarilll993l ). Understand- 
ing uj Cen may therefore assist our understanding of mass 
loss and chemical enrichment in other nucleated dSphs, as 
well as in the earlier, more metal-poor Universe. 



1.3 Individual stars 

Within this paper, we also present the first mid-infrared 
spectra of the two brightest, most IR-excessive stars in 
the cluster, V6 ( LEID 33062, ROA 162) a nd V42 (LEID 
44262, ROA 90). lGlass fc Feastl (| 19731 . 1 19771 ) measured the 
IR colours of V6, confirming it as a very bright IR source 
and showing it to be variable at near-IR wavelengths, 
with a possible L-band exc ess. Its period i s unc ertain, be- 
ing listed as 73.513 days in | Sawver Hoggl ( 1973) and 100- 
120 days in iDickens. Feast fc Lloyd Evans ( 19721 ). Classi- 
fied as an M4-5 emission line variable, it has a possi- 
ble radial velocity variation of up to 40 km s _1 with a 
mean vel ocity of +213.6±3.7 km s" 1 (Dickens et al. 1972; 
IWebbinkl [l981; quoted uncertainty ~13 km s" 1 ). Its tem- 
perature has been estimated at 3300-3600 K with log(q) 
of 0.0 (assuming M = 0.8 M B ; iPersson et al.|[l980l ; iFrogell 
1983). It contains relatively large amounts of H2O com- 
pared to the other cluster long- period variabl e s (LP Vs) 
and is CN and NH enhanced (|Cohen fc Belli Il986l ). It 
is also known to b e a TiO variable an d show s variable 
hydrogen emission (|Llovd Evans! Il983dl lal. [l986). V6 ap- 
pears to belong to t he metal-intermediate sub-population, 
with [Fe/Hl ~ -1.19 ilZinn fc Westjll984l; iNorris fc Da Costal 
ll995l; IVanture. Wallerstein fc Suntz cfi 2002). 

iFeastl (|196a ) implies V42 should be classified as a semi- 
regular variable of type SRd, but its spectrum is somewhat 
later than the F-K-type this implies. An Ml-2.5 emission 
line variable, it may also exhibit radial velocity variations 
between about +253 and +272 km s -1 (Dickens et al. 1972), 
though this variation may be largely attributable to insuffi- 
cient signal-to-noise. V42 may represent a star bridging the 
gap between SRd stars and emission-line LPVs. As an appar- 
ent fundamental-mode pulsator, its 148.64+0.03 day period 
(vL+00) is the longest in ui Cen and among the longest in 
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globular cluster variables as a group (|Clemendll997T l. Dick- 
ens et al. also report that the TiO bands weaken and hydro- 
gen emission lines are present only near o ptical photomet- 
ric maximum (see also lLlovd E vans 1983a, c), suggesting an 
additional opacity source. It is also CN and NH enhanced 
ijCohen fc Belli Il986h . ICacciari fc Freeman! (119831) estimate 
log(ff) = 0.5 and T = 3950 K. whereas iMenzies fc Whitelockl 
(| 1985h calculate a much lower temperature of 2818 K, based 
on IR data, and show V42 undergoes substantial variabil- 
ity - some ~60% of its mean luminosity - even in the L- 
band. No reliable information on the star's metallicity is 
available, although vL+07 place it at [Fe/H] ~ —1.25; a 
high-resolution optical spectrum shows significant H-a emis- 
sion, attributable to shocks p ropagating in the stellar wind 
IIMcDonald fc van Loon 2007. hereafter MvL07). 

The cluster also contains a number of other inter- 
esting objects. Five carbon stars have now been reported 
(vL+07). Several post-AGB stars are also known. Most fa- 
mously, Fehrenbach's Star (LEID 16018, HD 116745) ap- 
pears to have already undergone thorough mixing and mass 
loss dFehrenbach fc DuflotJ Il962l ; iDickens fc Powelll 1 19731 ; 
iGonzalez fc Walrersteinlll992T l. Another, VI (LEID 32029), 
is an irregularly-pulsating star with multiple periods that 
is thought to have undergone gas-dust separation, but 
not thought to have undergone third dredge-up (surface 
enrichment during t hermal pulses) and may thus be a 
post-early-AGB star dGonzalezl 1 19941 ; iMoehler et al.1 Il998l ; 
iThompson et al.ll200d 120071 ; vL+07). V29 (LEID 43105), 
V43 (LEID 39156), V48 (LEID 46162) and V92 (LEID 
26026) have also been suggested to be po st-AGB stars 
(|Gonzalezlli99i IGonzalez fc Warlersteinlll994T l. 

The remainder of the paper is organised as follows: Sec- 
tion 2 describes the SED input data and models; Section 3 
details the corrections to the photometry and the process of 
creation of the SEDs; Section 4 compares our observed tem- 
peratures and luminosities with those predicted by a variety 
of stellar evolution models, providing a new estimate of the 
distance and reddening to the cluster; Section 5 presents new 
mid-infrared spectra of V6 and V42, and discusses subse- 
quent estimation of mass loss from individual stars; Section 
6 discusses the implications of the dataset, including cal- 
culating the total mass-loss rate of the cluster, correlations 
with various stellar parameters, and the evolutionary nature 
of mass loss; finally, Section 7 presents our conclusions. 



2 THE INPUT DATASETS 

2.1 Literature photometry and variability data 

Our input list of objects was derived from the optical pho- 
tometry of vL+00. We have selected from this catalogue 
those stars that have at least one detection in the mid-IR 
Spitzer IRAC (3.6, 4.5, 5.8 and 8 /im) and MIPS (24 fim) 
imaging published in B+08. We have combined these data 
with photometry from 2MASS (|Skrutskie et al.ll2006l ). Of 
the resultant list of 6875 stars, 6018 are proper motion mem- 
bers, 1145 are radial velocity members in vL+07, and 1701 
have photometry at both 8 and 24 /jm. 
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Figure 1. Distribution of internal errors in temperature over the 
giant branch stars. Errors increase to typically ± ~1000 K for 
stars not covered by our hottest models (> 6500 K). 



2.2 The model spectra 

As many of our stars - particularly the cooler objects - de- 
part significantly from a blackbody spectrum, we have used 
model spectra to derive the stellar parameters. A grid of 
spectra were created using the Model Atm osphere in Radia- 
tive and Co nvective Scheme (marcs) code jGustafsso n et ail 
Il975ll2008f ) from 4000 K to 6500 K in steps of 250 K, with an 
additional dataset at 3500 K; gravities were sampled from 
log(g) = 0.0 to 3.5 in steps of 0.5. At each grid point, a 
further dimension in metallicity was sampled at [Z/H] = 
0.0 with solar abundances; at [Z/H] = -1.0 with [a/Fe] = 
+0.3; and at [Z/H] = -1.5 and -2.0, with [a/Fe] = +0.4. For 
each grid point, a synthetic spectrum was created with R = 
20 000 over the range 0.13 to 20 /im. Some models did not 
converge; these were replaced with a neighbouring conver- 
gent model that had the same temperature and the closest 
available combination of metallicity and gravity. 

The spectra were extrapolated beyond 20 /im using a 
Rayleigh- Jeans tail fit to the range 18-20 /im. This was 
used in prefer ence to the normally more-appropria te En- 
gelke function (|Engelkeill992l; iDecin fc Erikssonl feoO? ) as the 
Engelke function over-estimated the flux we observe by up 
to 20-40%, while the over-estimation from a Rayleigh- Jeans 
tail is only around 10-20%. Both a Rayleigh- Jeans tail or 
Engelke function are expected to tinder-estimate the flux 
we observe at 24 /im for our cooler stars, though the reverse 
may be true in warmer giants. The immediate reasons be- 
hind the apparent cwer-estimation these functions provide 
when compared to our data are not clear and we make an 
empirical correction in Section [5]3TT] 



3 SPECTRAL ENERGY DISTRIBUTIONS 

Firstly, the literature photometry and model spectra were 
converted into F v ( Janskys) . The models were then degraded 
in resolution by a factor of 100 to R = 200 to save on com- 
puting time and, using a cubic spline, interpolated in wave- 
length onto filter transmission data in the relevant filters. 
Comparing the photometric flux from the original spectrum 
and the interpolated spectrum, we find the differences to be 
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Figure 2. Distribution of internal errors in temperature over the 
cluster, containing only stars on the giant branches (here denned 
at T < 5500 K). Errors increase towards the cluster centre due to 
blending. The horizontal smearing is an interpolation artefact. 



■C 1% in all cases except _B-band, where the lower-resolution 
sampling of strong spectral lines lead to the systematic over- 
estimation of the flux by up to ~3% for the cooler stars. 
These are comparable to the errors present in the photome- 
try, and much less than the observed photometric variation 
of our mass-losing candidates, thus we do not consider it to 
be a major source of error in our analysis. An additional 
source of systematic error comes from uncertainties in the 
filter and imaging device responses, and instrumental pho- 
tometric zero points. Though difficult to quantify, we expect 
these effects to affect our photometry at the level of 3% or 
less and correspondingly systematically alter our tempera- 
ture and luminosities by ^$1%, as the errors will tend to 
average out among the different filters. 

A blackbody was fitted to the observed photometry to 
provide a first-order estimate of temperature and luminosity. 
From this we derive an initial value for \og(g), assuming a 
mass of 0.8 M . 

An initial grid of MARCS models was set up to coarsely 
define the temperature. For computing purposes, the grid 
was chosen to be between 4096 and 6144 K, in steps of 512 
(= 2 8 ) K. For stars falling outside this range (as determined 
by the fitting procedure), the grid was extended to 3072 K 
in the case of cooler stars, and to 18 442 K for warmer stars. 
The step size was chosen to provide the computationally- 
fastest determination of temperature, while ensuring the 
true best-fit temperature is reached. For each of these tem- 
peratures, a spectrum was produced from the closest two 
model grid temperatures by linearly interpolating in [Fe/H] 
and log(g). These two spectra (Fi and F2) were then com- 
bined into a single spectrum (-Fanal) by dividing the spec- 
tra by blackbodies of the corresponding grid point temper- 
atures, averaging the spectra, then multiplying by a black- 
body of the final required temperature. 

By comparing the interpolated spectrum with the orig- 
inal, we calculate that our interpolation method can accu- 
rately reproduce a model value for a photometric band at 
any temperature above 4000 K, with any gravity and metal- 
licity within our model grid to within 1%. For stars below 
4000 K, a lack of converging models at 3750 K and < 3500 



K means that the models will tend to over-estimate the pho- 
tometric flux at short wavelengths (B- and V-band) by per- 
haps 3-5% for the very coolest stars we sample, which is 
still (much) less than the expected stellar variability, differ- 
ential interstellar reddening or (in some cases) photometric 
uncertainty in this regime. 

Since ui Cen exhibits moderate but noticeable interstel- 
lar reddening, it was necessary to redden the model spec- 
trum appropriately. Estimates of the reddening towards the 
cluster vary considerably from E(B — V) = .09 to 0.15 mag 
|DiorgovskillT993l : lHarrislll996l ; |Kovacsl l2002). A differential 
reddening across th e cluster of AE(B — V ) ^ 0.04 mag has 
also been inferred (|Norris fc Besselll fl975l ; ICalamida et al.l 
120051 ; vL+07). For our analysis, we adopted a reddening of 
E(B — V) — 0.08 mag across the entire cluster to better 
fit the evolutionary is ochrones (see Section 2J. Assuming 
A V /E(B - V) « 3.05 (|WhittetJll992T ), this yields an extinc- 
tion of Ay ~ 0.24 mag. 

To compute the reddening as a functio n of wa velength, 
we have used the relationship from iDraind l|l989l ). with an 
empirical correction for wavelengths of < 1 /im based on 
Draine's Figu re 1. Using E(J - K ) = 0.52E(B - V) (after 
IWhittet|[l992 ), this becomes: 



E{B~V) 



1.248A- 1 ' 75 



for A ^ 1 /im 
for A ^ 1 /im 



■(1) 



By examining the effective temperature and luminosity 
at different reddening factors on a subset of stars, we can 
deduce the size of this correction. We find that the corre- 
sponding temperature and luminosity, T' and L' (in Kelvin 
and solar luminosities), are well approximated (for an arbi- 
trary E(B — V), denoted e) by: 



T — T + 



0.08 
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0.08 
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13200 



-0.155 



1.5300 



(2) 



(3) 



where L and T are the luminosity and temperature for 
E{B — V) = 0.08 mag and d is the distance in parsecs. 

Draine estimates a 10% error in his relationship, which 
would alter our corresponding corrections by a similar 
amoun t. Our composit e relation described above differs from 
that of lWhittetl (|l992T ) by 10 - 15% over most of the spec- 
trum, suggesting a similar inherent uncertainty in our red- 
dening correction. The uncertainty in E(B — V) affects the 
shorter wavelengths the most, and is typically greater than 
that produced by our filter convolution and interpolation 
above. Differential reddening across the cluster thus leads to 
a scatter of up to ±50 — 100 K in temperature and ±1 — 4% 
in luminosity. Uncertainty in the average extinction to the 
cluster, which is of the same order, is the largest systematic 
error in our temperature and luminosity estimates. 

We do not consider circumstellar reddening here, as it 
only becomes significant in individual stars for the most ex- 
treme cases, where the effect is dwarfed (by a factor of ~10) 
by the variability due to pulsation at any particular wave- 
length. This will not significantly affect our determination 
of distance or extinction to the cluster, but may lead to a 
slight under-estimation of the temperature and luminosity 
of the most enshrouded stars. 

Once appropriately reddened, the spectrum was con- 
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volved with the filter responses to obtain the expected flux 
density received in that filter. A true \ 2 statistic was im- 
possible due to incomplete error information; a pseudo-x 2 
statistic was calculated between the observed and expected 
flux densities, setting all observed photometric errors to have 
equal weight. The 24-/^m data were not included in this cal- 
culation as we aim to fit the underlying photosphere, not 
the photosphere plus the wind's dust component. 

The process was repeated for each of the points in 
the one- dimensional temperature grid and the x minimum 
identified. The value of log(gr) was re-determined for the 
temperature of the x 2 minimum (T, 
determined at half the grid spacing, T, 
X 2 minimum re-defined. This process was iteratively carried 
out to refine the temperature to a precision of 1 K. The re- 
sulting x 2 distribution suggests that the internal accuracy 
is typically ±70 K, but increases towards lower luminosi- 
ties and higher temperatures. Note that while colour- T en 
transformations often give more precise values for T e ff, our 
procedure yields the most consistent values for T c s and lu- 
minosity, taking into account the full observed SED. We plot 
the temperature errors in Fig. [1] 

In Table [1] we take the average and standard deviation 
of the difference between the MARCS modelled flux and the 
observed flux in each photometric band. We have limited 
this to stars with T e g < 6500 K, due to our lack of models 
at higher temperatures. Note also that the accuracy of our 
temperatures is mostly determined by the accuracy of our 
B and V-band flux. 

We include an interpolated map of the temperature er- 
rors in Fig. [2 showing they increase toward the cluster cen- 
tre. The giant-branch stars with larger temperature errors 
(<;200 K) are caused by errant photometry, particularly in 
the B- and V-bands, due to source blending in the dense 
cluster core. As illustrated in Fig. [1] there is only a very 
weak dependence between the magnitude of the error and 
the stellar temperature (and thus luminosity if the star is 
on the giant branch). The frequency of errors due to blend- 
ing does increase as one heads to lower temperatures and 
luminosities. Crucially, however, the very presence of these 
errors allows us to determine which stars suffer from blend- 
ing. We find that, on the upper giant branch, where mass 
loss is taking place, the only stars to suffer from substantial 
errors (130-330 K) are known variables (specifically V42, 
V6, V152 and V148), whose temperatures are expected to 
be uncertain due to their inherent variability. 

Literature distance e stimates vary from 4. 8 kpc 
(vL+00) to 5.52±0.13 kpc (|Del Principe et all [20061' ) with 
5.3 kpc being the m edian I Peterson! 19931 ; Harri s 1996; 
iThompson et al.l200 if ). We adopted 5.0 kpc as an initial esti- 
mate, based on evolutionary isochrones (Section 3J). Assum- 
ing this distance, a luminosity was calculated by integrating 
the final model spectrum. Due to the limits of integration 
(>130 nm) determined by the model spectra, we note that 
the luminosities of stars with considerable flux at A < 130 
nm (i.e. stars with T 2> 8000 K) will likely be more lumi- 
nous than we have listed here, though this is not certain as 
our stellar models become progressively less reliable beyond 
6500 K. The final stellar parameters are listed in Table [2] 



Table 1. Average differences between modelled and observed 
fluxes, and standard deviations of those differences, for stars be- 
low 6500 K. 



Band 


Number 


Average 


Standard 




of stars 


difference 


deviation 


B 


5867 


+5.1% 


9.9% 


V 


5861 


-6.9% 


9.9% 


J 


5645 


+0.9% 


7.2% 


H 


5598 


-2.2% 


6.7% 


K 


5585 


+4.2% 


8.8% 


3.6/jm 


5463 


+2.6% 


8.4% 


4.5/^m 


5673 


-0.6% 


8.3% 


5.8/im 


5319 


+0.4% 


14.5% 


8.0/rni 


5759 


-1.3% 


12.6% 



Table 2. Derived stellar parameters (full table on-line) ordered 
in decreasing luminosity, assuming E(B — V) = 0.08 mag and a 
distance of 5.0 kpc, including Leiden Identifier and proper mo- 
tion (PM) percentage likelihood of membership from vL+00, and 
radial velocity from vL+07, where available. 
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A physical Hertzsprung-Russell diagram (HRD) is 
shown in Fig. [3] The RGB and AGB can clearly be distin- 
guished up to ~150 L.0, with the Horizontal Branch (HB) 
extending towards high temperatures at ~50 Lq. The ex- 
pected errors suggest the majority of the spread in the di- 
agram is real. The cooler giant branches should therefore 
represent the more metal-rich objects. Post-AGB stars are 
present at higher luminosities towards the warm side of the 
giant branch. The RGB-a stars are visible to the right of the 
main RGB. Interestingly, they do not appear to extend to 
the RGB tip, but as they are so few in number, this could 
merely be a stochastic effect. 



4 COMPARISONS WITH STELLAR 
ISOCHRONES 

With an HRD of this quality, a comparison with stellar evo- 
lution models can yield accurate determinations of the clus- 
ter's parameters, such as distance and reddening (we assume 
a fixed age of «12 Gyr as this is better estimated using the 
main-sequence turnoff). These parameters are important for 
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Figure 3. A physical HRD for the giant branch of uj Cen. Dots show proper motion members (>50% probability); squares show objects 
that are (also) confirmed radial velocity members (ulsr > 100 km s — 1 ). Representative statistical errors for individual object s are s hown 
in the lower- and upper-right corners. The limits to the RR Lyrae instability strip are taken from lBono, Caputo fc Marconil (jl995). 



both calibrating our luminosities and estimating mass-loss 
rates from individual stars. Evolutionary tracks and, in par- 
ticular, the rate of stellar evolution, can already yield some 
constraint on mass loss on the RGB. 



4.1 Padova isochrones 

Fig. [4] shows Padova stellar isochrones from iMarigo et alj 
(2008). The isochron es are for zero extinction, c o ntain ing the 
dust solutions from iBressan. Granato fc Silval (119981), wit h 
a log-normal initial mass function from IChabrierTl|200ll) . 
They are shown for the metallicities and age s suggested by 
lOriglia et all (120031 ) and iHilker et al.1 (|2004h . namely: 12.1 
Gyr at [Fe/H] = -1.6 for the metal-rich population, 10.5 
Gyr at [Fe/H] = -1.2 for the metal-intermediate population 
and 9.4 Gyr at [Fe/H] = -0.7 for the metal-rich population. 
Using these isochrones, we fitted the extinction and distance 
as E(B -V) = 0.09 ± 0.01 mag and 4800 ± 150 pc, respec- 
tively. 

We stress that these two parameters are somewhat cor- 
related, and the errors given are based on single-parameter 
errors only. In this case, the distance is fitted primarily by 
the location of the zero-age horizontal branch (ZAHB), and 
the isochrones are then matched to the data by the correct 



reddening. While we have not undertaken a full quantitative 
analysis, comparing the two panels in Fig.|4j we can see that 
a fit with a distance of 5000 pc and E(B — V) of 0.08 mag 
also falls within visually plausible errors, however this does 
not reproduce the early- AGB and early-RGB so well. 

Despite efforts to fit the distance and reddening, the 
Padova isochrones still do not provide a good fit to the slope 
of the RGB or AGB. We suggest that this is inherent in the 
models for three reasons: 

• on the basis of our above analysis, circumstellar dust 
starts reddening the stars at lower luminosities than their 
models predict, at around 1000 L© (c.f. Fig.[T7J - this occurs 
both on the RGB and AGB; 

• RGB mass loss is not included in the Padova models 
(the accumulated mass loss is assumed to happen 'instan- 
taneously' at the RGB tip): stars which have lost mass will 
have larger radii and thus be cooler than the models predict; 

• the [cv/Fe] enhancement in the models is incorrect 
(Marigo et al. currently do not provide a mechanism with 
which to alter this) - increased [a/Fe] will make the stars 
cooler than the models. 

It also appears that the mass-loss rates from IMarigo et al.l 
(2008) for the RGB are not high enough for a large number 
of stars, as they do not reproduce the blueward extent of the 
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Figu re 4. HRD of cluster members, including isochrones from 
iMaiigo et all ll2008f) . The top panel shows isochrones for metal- 
licities of [Fe/H] = -1.6, -1.2 and -0.7 (from left to right), as- 
suming a distance of 5000 pc. The arrows show the approximate 
reddening correction (taking E(B — V) = 0.08 mag) for that tem- 
perature. The bottom panel shows the same diagram, but with 
an approximate correction to E(B — V) = 0.09 mag and a dis- 
tance of 4800 pc. The isochrones are described in the text, with 
metallicity increasing from left to right. 



HB that results from low mantle masses. We can therefore 
assume that, based on the Padova models, a star with an 
initial mass of 0.843 Mq, must typically lose >0.117 Mq on 
the RGB. 



4.2 Dartmouth isochrones and ZAHB models 

Fig.[5]shows a si milar plot, but for iso chrones from the Dart- 
mouth database (|Dotter et al.|[20oii ). In the top panel, we 
take a 12.1 Gyr model at [Fe/H] = -1.62, a standard he- 
lium abundance of Y — 0.245 + 1.5Z and an [a/Fe] ratio 
scaled to solar abundances. We also include a HB model at 
M = 0.70 M Q , implying again that ~0.11 M Q is lost on the 
giant branch (the RGB tip corresponds to an initial stellar 
mass of 0.81 Mq in this model). 

By applying a reddening correction to E(B — V) = 
0.08 ± 0.01 mag and a distance of 4900 ± 100 pc, we can 
yield a good match to the isochrones: the RGB is matched 
nearly exactly, as is the AGB. The HB morphology is not 
well fit by solar-scaled helium abundance and [a/Fe]. In the 
middle panel of Fig. [5] we show the effect of varying mantle 
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Figure 5. As Fig. [4] using RGB isochr ones and HB model s gen- 
crated using the Dartmouth database jPotter et al]|2008h . The 
top panel shows the same three RGB isochrones as Fig. [4] with a 
solar-abundance-scaled HB track at 0.70 Mq. The middle panel 
shows the same diagram, using only the 12.1 Gyr isochrone, but 
with an approximate correction to E(B — V) = 0.08 mag and a 
distance of 4900 pc, with HB models at (top to bottom) 0.65, 0.70, 
0.75 and 0.81 Mq. The bottom panel is as the middle panel, but 
with [a/Fe] enhancements of (bottom to top) 0.0, +0.2, +0.4 and 
+0.6; the RGB isochrones are shown for [a/Fe] = 0.0 (left) and 
[a/Fe] = +0.4 (right), both have helium abundances of Y = 0.33. 
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Figure 6. As Fig. [5] middle panel, showing a Dartmouth model 
covering the early AGB for a stellar mass of 0.65 Mq, at [Fe/H] 
= -1.62 and [ct/Fe] = +0.2, with reddening (from left to right) 
of E(B - V) = 0.11, 0.10, 0.09, 0.08 and 0.07 mag. The data arc 
shown at E(B - V) = 0.08 mag and d = 4900 pc. 

mass at solar [a/Fe] with Y — 0.33. In the bottom panel, 
we show the effect of fixed stellar mass with varying [a/Fe]: 
this affects the mass of the stellar core which, for fixed stel- 
lar mass, has the reverse effect on the mantle mass. We note 
that the errors given above are in that parameter only, how- 
ever we can rule out a fit at E(B — V) = 0.07 mag and 
d = 5100 pc, or at E(B - V) = 0.11 mag and d = 4500 pc 
by examining the locations of the start of the AGB and the 
tip of the RGB. 

It seems clear that variation of the mantle mass can 
account for the observed spread in the HB of u> Cen. This 
arises from a combination of varying core mass due to intrin- 
sic metallicity and abundance differences. It can also arise 
from different efficiencies of mass loss on the RGB, which 
could well be due to the original metallicity and abundance 
differences themselves. At first sight, it is not clear whether 
varying mass loss or varying core mass is the primary fac- 
tor. However, the spread of the early AGB is very narrow 
(~10-15% in luminosity at a fixed temperature, 50-90 K in 
temperature at a fixed luminosity), most of which is due to 
temperature and luminosity errors of order 50-70 K, with 
small contributions from differential reddening across the 
cluster (implying, from Eq. J2|, AE(B - V) < 0.02 mag) 
and varying distance within the cluster. It would therefore 
seem that we need a range of models that produce almost 
zero spread in the early AGB, a condition satisfied by chang- 
ing the stellar mass, rather than [a/Fe]. 

We can further estimate the maximum value of AE(B — 
V) across the cluster. Fig. [6] shows the effect of different 
values of E(B — V) on the early AGB. In the absence 
of other variations leading to the dispersion of the AGB 
(which include photometric errors, and [a/Fe], metallicity 
and distance variation), we can assume that the variation of 
E(B - V) across the cluster, AE(B - V), is < 0.04 mag. 

By examining the HB (Fig. [7]), we can approximate the 
mass loss on the RGB. It would appear that the majority 
of the HB stars lie between the 0.60 and 0.65 Mq models 
in both temperature and luminosity, but with some more 
closely following the higher-mass models. We must stress 
that our models do not cover the HB: while we expect tem- 
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Figure 7. As Fig. \E\ middle panel, showing the HB for stellar 
masses of (top to bottom) 0.60, 0.65, 0.70 and 0.81 M , at [Fe/H] 
= -1.62 and [a/Fe] = +0.4. 

peratures to be broadly correct (within a few hundred to 
a thousand Kelvin), luminosities may be too low. We may 
also not cover the blue extent of the HB in its entirety, as 
bluer HB stars will be fainter at Spitzer's IR wavelengths 
than their redder counterparts. 

A stellar mass of around 0.62-0.64 Mq suggests that 
0.18-0.20 Mq is lost on the RGB, which is in keeping with 
iDotter I (|2008l ). who suggests that the total mass lost on the 
RGB does not greatly depend on metallicity above [Fe/H] 
~ -2. 

Globular cluster white dwarf masses (e.g. lMoehler et ail 

|2004|) suggest that ~0.3 Mq must be lost on the RGB and 
AGB combined (the core mass in this case is 0.49 Mq), 
implying that around two thirds of the mass loss occurs on 
the RGB in ui Cen. Interestingly, we also expect about a 
third of the mass loss from the cluster to be in the form of 
chromospheric winds (see Section I6.5.2|l , implying at least 
half the dusty mass loss occurs on the RGB. 

Some concerns have been raised that the HB models 
from the Dartmouth models may be marginally too warm 
for low- metallicity stars (A. Dotter, private communication). 
Nevertheless, given the much more precise fit to the data, 
we would ad vocate the Da r tmou th model values over those 
derived from lMarigo et al.l (|2008l ). 



4.3 Victoria-Regina isochrones and ZAHB models 

We have performed a similar analysis using the 
Victoria-Regina (VR) isochrone and Z AHB models of 
IVandenBerg. Bergbusch fc Dowierl l|2006h . which we show 
in Fig. Here, we take the models at [a/Fe] = +0.0, +0.2, 
+0.4 for [Fe/H] = -1.61 at 12 Gyr. No values of reddening 
or distance, nor any sensible value of age can accurately 
match the VR isochrone on the RGB to our data for any 
value of [a/Fe]. In order to achieve a fit, we must either 
increase the temperature of the VR model, or decrease the 
temperatures of our stellar data, by around 0.01 dex (2.3%) 
in the case of [a/Fe] ~ or 0.02 dex (4.7%) in the case 
of [a/Fe] ~ +0.4. The value of [a/Fe] is suggest ed to be 
aroun d +0.3 for cluster stars with [Fe/H] —1 ()Pancinol 
2003). Under this assumption, we find a reddening of 
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Figure 8. As Fig. [4] using Victoria- Regina isoch r ones a nd ZAHB 
models from IVandenBerg, Bergbusch & Dowlcr (200(|). The top 
panel shows a 12 Gyr isochrone and ZAHB model for [Fe/H] = — 
1.61, with [a/Fe] = +0.0 (solid line), +0.2 (dashed line) and +0.4 
(dotted line). The bottom panel shows the fitted model at [a/Fe] 
= +0.4, still at 12 Gyr, with E(B-V) = 0.09 mag and d = 4800 
pc, with the model temperatures increased by 4.7% (0.02 dex). 



around E(B -V) = 0.08 mag and a distance of 4800 ± 100 
pc. 

On the assumption that the ZAHB location for the ma- 
jority of metal-poor stars is well-represented by the clump 
around 8500-11 000 K, and that [a/Fe] « +0.3, we find that 
most stars have a ZAHB mass of around 0.61+0.02 Mq (Fig. 
[9]). Obviously, this issue is clouded by HB evolution and lack 
of fully accurate measurements of temperatures and lumi- 
nosities in this area, but the figure is broadly consistent with 
that from the Dartmouth models above, predicting only a 
slightly lower stellar mass. The 12 Gyr VR model predicts 
that the initial mass of RGB tip stars is 0.85 Mq and that 
the core mass is 0.49 Mq at a — +0.3. This suggests that 
«0.24 M is lost on the RGB and ~0.12 M on the AGB, if 
the entire envelope is ejected. This ratio of 2:1 is consistent 
with that derived from the Dartmouth models. 



4.4 BaSTI isochrones 

We show a fit using th e BaSTI models ('A ba g of stellar 
tracks and isochrones'; l)Pietrinferni et al1l2006h ) in Fig. \W\ 
(top panel). As before, we take a 12 Gyr model at [Fe/H] 
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Figure 9. As Fig.0 with d = 4750 pc and E(B-V) = 0.08 mag, 
showing the ZAHB models from the Victoria-Regina models, at 
[Fe/H] = -1.61 and [a/Fe] = +0.0, +0.2 and +0.4 (as in Fig. 
[8)l . Thicker lines show the stellar masses on the ZAHB at 0.575, 
0.600, 0.625 and 0.650 Mq. 
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Figure 10. As Fig. [4] using BaSTI isochrones. The top panel 
shows the same three RGB isochrones as Fig. [4] with the en- 
enhanced AGB-extended BaSTI model for Z = 0.001 ([Fe/H] 
= —1.62), at 12 Gyr. The bottom panel shows the fitted model, 
at 9.5 Gyr, with E(B-V) = 0.06 mag and d = 4900 pc. The solid 
lines assume Reimers' r] = 0.4; the dotted lines show the AGB for 
f] = 0.2. 
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= -1.62 at a helium content ol Y = 0.246. We u se the a- 
element enhanced models with r\ = 0.2 and 0.4 |Reimersl 
Il975l ). inco rporating BaST I 's syn thetic AGB evolutionary 
tracks from llben fc Truranl (ll978T ). 

We note that, although the synthetic HB and AGB data 
fit the empirical cluster data well, the RGB is significantly 
cooler and under- luminous in the BaSTI model. The only 
way we can find to rectify the large gap between the early 
AGB and RGB is by assuming the cluster is much younger 
than has previously been calculated, at an age of < 9.5 Gyr. 
Taking models at this age, we found the reddening to only 
E(B - V) = 0.06 mag and d = 4900 pc. 

The young age of the fitted model is at odds not only 
with the luminosities of the RGB and AGB tips, but with 
previous age estimations of typically 11-14 Gyr, derived 
using a variety of data and methods including isochrone 
fitting to colou r-magnitude diagrams including stars down 



to ~ 0.15 Mm (iThompson et alj|200ll; iKaluznv et al 
Chabover fc KrausslbOOsi r iPancin o 2003; Platais et al 



2002 



2003 



Hilker et al.ll2004 ). BaSTI models this old do not have the 



narrow separation between the RGB and AGB. Another im- 
plication of this is that mass-loss efficiency must range from 
77 w 0.2 to n » 0.4 (probably w 0.7 to 0.8 and perhaps 
higher for the Extreme HB, which we do not model). This 
level of mass-loss efficiency is not considered to be impossi- 
ble, and ver y similar values hav e been suggested for ui Cen 
in the past (|D'Cruz et al.lll995T ). 



4.5 Summary of isochrone fitting 

In summary, we find that the Dartmouth mode ls provide the 
best fi t to our data, although the models of iMarigo et al.l 
(2008) may provide a similarly good fit if the mass-loss ef- 
ficiency could be increased. The BaSTI models provide too 
small a temperature difference between the early RGB and 
AGB, which would appear unlikely to be due to a systematic 
error in our data, and can consequently only fit our data with 
a significantly decreased age. The Victoria- Regina models 
have an effective temperature which clearly differs from our 
data by up to 5% in both the isochrones and the ZAHB 
models. 

Taken in combination, the models consistently suggest a 
distance to the cluster of around 4850±200 pc, and an aver- 
age reddening to the cluster of E(B — V) = 0.08 ±0.02 mag. 
These values are both at the lower limit of those derived in 
the literature, which roughly span the ranges of 5160±360 
pc and E(B - V) = 0.12 ± 0.03 mag (see Section |3] for ref- 
erences). The distance is very close to the 4740±160 pc in 
vL+00, whose data we use here. However, we note that this 
was derived using velocity dispersion data, not the photom- 
etry we use here, and is thus nearly independent. We cannot 
entirely rule out a distance of 5000 pc, though this fits the 
start of the early AGB less well, and requires a reddening 
of only E(B — V) = 0.05 mag. Conversely, a reddening of 
E{B — V) — 0.12 mag would suggest a distance of ~4600 
pc. 

The potential for systematic error in our temperatures 
and luminosities (Section [3} is relatively small for stars 
within the range 4000-6500 K, and increases as one moves 
away from this region. We estimate that, excluding uncer- 
tainties in distance, and average and differential reddening, 
systematic differences on the RGB and AGB can be lim- 



ited in temperature to < 1% and in luminosity to i; 5%. 
Very red stars and HB stars have the potential for higher 
systematic errors due to the lack of MARCS models covering 
these temperature ranges. This has the potential to make 
a systematic offset to our distance of up to ±120 pc and 
of up to ±0.015 mag in E(B - V). The 10-15% discrep- 
ancy in our reddening correction found in Section [3] and the 
assumption that the foreground ISM follows the approxima- 
tion Av/E(B — V) ~ 3.05 probably inflates the systematic 
error in E(B — V) to 0.020 mag. Treatment of convection 
by the stellar evolution models may also cause a systematic 
offset in the distance and reddening we determine to the 
cluster. However, all four models include physics describing 
convective core overshoot. The exact physics employed does 
not appear to have made a significiant difference to our re- 
sults, as the distance we estimate to the cluster is roughly 
the same with all four models. 

As interstellar extinction has a larger effect at bluer 
wavelengths, the alteration of E(B — V) does not much affect 
the values we derive for our mass-losing stars, which are 
typically cooler than the giant branch, although a distance 
of 4850 pc would mean that their luminosities would still 
decrease by 6% due to the distance correction. The net effect 
would be a decrease in our total mass-loss rates by 4.5%, 
though this is included in the errors we list in Sections 15.21 
fcl6~2l 

From arguments stemming from the location of evolu- 
tionary sequences on the HRD, we conclude that most of 
a cluster star's mass loss must occur near the tip of the 
RGB, rather than on the AGB. The HB morphology of the 
cluster requires that typically 0.20-0.25 Mq of stellar at- 
mosphere is lost from the (initially) 0.85 Mq stars on the 
RGB. Stars then attempt to lose a further 0.05 Mq or more 
on the AGB, but many may have insufficient atmospheric 
mass to lose and thus become post-early- AGB st ars or even 
AGB-manque stars ('failed' AGB stars - see, e.g. lO'Connelll 
Il999l , Section 6.2). That mass loss can work this efficiently at 
metallicities as low as [Fe/H] ~ —1.6 suggests that metal- 
licity and mass loss have only a weak dependence in this 
metallicity regime. 



5 DERIVING MASS-LOSS RATES 

5.1 Mid-IR spectra and additional data of V6 and 
V42 

Spectra were taken of the two most extreme red giants in 
the cluster: V6 (LEID 33062, ROA 162) and V42 (LEID 
44262, ROA 90). Taken wit h the T-ReCS spectrograp h on 
the Gemini South Telescope (|De Buizer fc FisheifeoOST ) ■ the 
data are between 8 and 13 fj,m and have a resolving power 
of R — 300 — 600, along with iV'-band (11 /im) acquisition 
images, on each of the nights of 2007 August 16 & 18, for 
V6; and on 2007 August 27 & 28, for V42. Two spectra of 
2.5 hours integration were taken each night, using the low 
resolution 10-pim grating, along with the 0.72" slit. The full- 
width half-maximum image width of the standard star 6 Cen 
in the iV'-band acquisition images was ~0.44". 

The spectra were reduced using the iraf packag 



1 IRAF is distributed by the National Optical Astronomy Obser- 
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Figure 11. Flux-density-calibrated Gemini T-ReCS spectra of 
the two extreme red giant stars lu Cen V6 and V42, overlaid with 
a running boxcar average of 10 pixels. 



signed for T-ReCS using the default settings. No wavelength 
offsets were observed, and differences among the spectra 
were small, so a flux-density-weighted average was taken 
to provide a single, high-signal-to-noise spectrum for each 
target. An absolute flux density for each target was esti- 
mated by comparing the acquisition images of the target 
stars and 6 Cen using aperture photometry. Assuming an 
iV'-band flux density of 45.4 Jy for 8 CerQ a flux density 
of 203±12 mjy was calculated for V6 and 91±16 mjy for 
V42. The acquisition field for V42 also included the mid-M- 
type irregularly- variable red giant star V152 (LEID 44277), 
which is a confirmed cluster member (vL+00, vL+07). By 
comparison to V42, we estimate its iV'-band flux density is 
36±8 mjy. Comparison with the 8- and 24-/jm Spitzer data 
shows this value to be consistent with a mid-infrared ex- 
cess, but the error on the Gemini photometry is too large 
to constrain whether a 10-fim silicate feature is present in 
V42. 

The flux-density-calibrated spectra for both stars are 
shown in Fig. 1111 Clear fringing in the spectrum of V6 can 
be seen. Attempts to correct for this by moving each spec- 
trum by an integer number of pixels failed to improve this. 
We note that many of the 'features' in the spectra are also 
the result of imperfect atmospheric correction, though the 
presence of 9.5-pm emission in V6 is certain. 

In the case of V6, additional o p tical p h otometry was 



in the case or Vb, a dditional op tical p h otometry was 
available from Cannon fc Stobid (| 19731 ) . iLlovd Evansl 
dl983bh. and IClementl i| 19971 ). Near-IR data from 
iPersson et ail (|l980l ) were also used. Of particular use 
were the wo rks by Dickens et al. (1972), and ldass fc Feasti 
l|l973l . Il977l ). which not only provide flux densities, but 
variability amplitudes, allowing us to see the temperature 
and luminosity changes in the star from U- to L-band. We 
have previously obtained a 2dF optical spectrum for this 
star (published in vL+07), where it was determined to have 
a temperature of ~3500 K (the lowest-temperature model). 




Wavelength (microns) 



Figure 12. SED models of V6 (top) and V42 (bottom), show- 
ing MARCS models smoothed to R = 200 (solid lines), smoothed 
Gemini T-ReCS spectra (dotted lines) and literature photometry 
(crosses with error bars). The error bars in the latter reflect either 
the photometric variability or the measurement error (see text for 
details). The dashed curve shows the corresponding DUSTY model, 
fitted to the Gemini spectrum and 24- /im Spitzer datum in the 
case of V6 and the literature photometry, including Spitzer data, 
in the case of V42. 



For V42, the same papers provided information on flux 
densities and variability. Additional data were available from 
the DEN IS survey (Deep Ne ar Infrared Survey of the South- 
ern Sky; iFouque et al1l2000l ); iMenzies fc Whitelockl i|l985l ). 
which includes variability inf ormation; and 9.6-, 10- and 12- 
[im photometric data po ints l|Origlia. Ferraro. fc Fusi Peccil 
Il995l ; lOriglia et al.|[2002l , hereafter O+02). In addition to a 
2dF spectrum, we also have an optical VLT/UVES spec- 
trum for V42 (published in MvL07), though our previous 
temperature estimates have not been particularly reliable 
as a result of veiling by emission dissipated in the pulsation 
shock. The resulting SEDs for V6 and V42 are shown in Fig. 

rn 



vatorics (NOAO), which is operated by the Association of Univer- 
sities for Research in Astronomy, Inc., under co-operative agree- 
ment with the National Science Foundation. 

2 Extrapolated from http:/ /www. gemini.edu/sciops/instruments/ 
miri/T-ReCSBrightStan.txt 



5.2 Deriving a mass-loss rate for V6 and V42 

5.2.1 Mass loss of V6 from Gemini spectroscopy 

Estimates of the mass-loss rate for the two stars V6 
and V42 were determined using the dusty modelling 
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code l|Nenkova, Ivezic fc Elitzuilll999h . Our Gemini T-ReCS 
spectrum of V6 (Fig. shows clear emission from sili- 
cate dust grains around 10 /im, typical of well-developed 
dusty winds seen aroun d AGB Mira vari ables in the So- 
lar Neighbourhood (e.g. ISpeck et alJl200Ch . The spectrum 
is not well fit by th e classic 'astronomical' silicates alone 
( Dra ine fc Lee] 1 1984) , especially the 9.5-fim peak, regard- 
less of the dust properties used. Our best fi t was for 65% 
(by n umber of grains) 'astronomical' silicates (I Draine fc Led 
11984 ). 15% compact AI2O3 (optical constants from Jena 
databasjfl), and 1 0% each of glassy and crystalline silicates 
l|Jager et al.ll 1994 ) (Fig. [HI). Although we have used these to 
produce the best fit, we do not claim that these are necessar- 
ily the particular species of silicates involved, nor that they 
are present in these proportions. We assume a radiatively- 
driven wind and a standard MRN grain size distribution 
l|Mathis. Rumpl fc Nordsieckl [l977l) . where the number of 
grains of size a is given by n(a) = a~ q , where we here take 
q — 3.5 over the range a = 0.005 — 0.05 /im. We derive an 
inner edge to the dust envelope at a temperature of 650 K, 
and a V-band optical depth of Ty = 0.08 for V6. This dust 
temperature is relatively cool, and could potentially be re- 
lated to mid-IR variability, which we discuss in Section r5.2.2l 
An increase in the dust temperature to a more typical 1000 
K would decrease the mass-loss rate by a factor of nearly 
two (Section [5321 . 

The maximum grain size also has a weak effect on 
the dust mass-loss rate. Here, we have chosen 0.05 /um, 
which we believe is representative of this metal-poor en- 
vironment, which may be subject to substantial UV flux 
from other cluster members. This also provides a marginally 
bett er fit to the region around 11 /i m in our V6 spectrum 
(c.f. IVoshchinnikov fc Hennind I2008I . Fig. 1). Our value is, 
however, less than the usual value of 0.1 ^m, or higher , 
assumed in the literature (e.g. iPapoular fc Pegourie]|l983h . 
Increasing the maximum grain size to 0.1 /im would de- 
crease the mass- loss rate by ~20%, a factor that does not 
change much as one further increases the parameter. While 
theoretically possible to measure, it is very difficult to con- 
strain the grain sizes involved purely f rom the SED alone 
|Carciofi. Biorkman fc Magalhaesll2004 ). 

The dust mass-loss rate and total mass-loss rate from a 
star a re related to the dusty output b y the following equa- 
tions l|Nenkova. Ivezic fc Elitzurlll999l ): 



M dus 



Mdusty / _L_\ 3/i (JL\ 
200 U0 4 / V200/ 



-1/2 



(f) 



1/2 



M = M dust ip, 



(4) 
(5) 



where L is the luminosity in solar units, ip is the gas-to-dust 
ratio and p B is the bulk grain densi ty in g cm " 3 (we assume 
it to be 3 g cm" 3 for silicates; e.g. lSuhlll999l ). 

Eq. (0 holds when condensation into grains occurs effi- 
ciently. We expect much of the condensible material to form 
silicate grains, as iron will typically form inclusions, rather 
than competing with silicon. We therefore expect silicate 
grain growth to scale inversely linearly with the metallic- 
ity, and specifically the abundance of silicon. Eq. ([5} may 
not extend to stars where chromospheric mass-loss is im- 
portant and dust production is not efficient, where the dust- 

3 http:/ /www. astro. uni-jena.de/Users/database/entry.html 



to-gas ratio is expected to follow a square-law dependence 
on metallicity (jvan Loonll2000l . [2006l) . We assume here that 
ij) = 200 for solar c omposition stars, which for V6 we take 
as [F e/H] = -1.19 (|Zinn fc Westlll984 ; iNorris fc Da Costal 
Il995l ; Vanture et al. 2002). 

These mass-loss rates are uncertain due to: 

• 30% uncertainty in the dust mass-loss rate from internal 
errors in DUSTY; 

• J; 30% uncertainty in the dust mass-loss rate due to 
uncertain dust chemistry, temperature and grain size; 

• 11% and 6% uncertainty in the photometric excess 
above the SED at 8 and 24 pm, respectively; 

• ~ 3% uncertainty in the dust mass-loss rate due to in- 
accuracies in the stellar temperature; 

• J; 5% uncertainty due to error in the integrated flux 
density from the SED; 

• 22% uncertainty due to error in luminosity; 

• ~ 5% uncertainty due to error in the reddening; 

• ~ 12% uncertainty due to error in the distance. 

When added in quadrature, this gives a total error of ^ 
52% in the dust mass-loss rate from V6, and is also based on 
the following assumptions, which we will discuss in Sections 
[532l fc[Q1 

• that the mid-IR spectrum V6 has not varied substan- 
tially between the 8- and 24-/im exposures; 

• that V oc 10" [Fc/H1 ; 

• that the wind velocity scales as prescribed in the dusty 
code: v oc L 1/4 V>" 1/2 ; 

• that the dust is coupled to the gas. 

We thus find that for V6, under the above assumptions, 
Mdust = 3.8 ±?;g xlO" 10 M Q yr" 1 and M total = 1.1 ±g;| 
xlO" 6 M yr" 1 . 

This total mass-loss rate for V6 is comparable to 
our earlier measurements in 47 Tuc (|van Loon et al.l 
2006), though it is generally higher than more metal- 
rich, solar-mass Mira vari ables in the Solar Neighbourhood 
|jura fc Kleinm ann 1992). It is also slightly inflated from 
those derived in B+08 purely on the basis of the Spitzer 
data and the models of iGroenewegenl (|2006l ). We note that 
the above values and quoted errors depend on the invariance 
of V6 between the acquisition times of the Spitzer IRAC and 
MIPS photometry, and of the Gemini T-ReCS spectrum. 



5.2.2 Mass loss of V 42 from Gemini spectroscopy 

The mass-loss rate of V42 is more complicated to derive, as 
there are considerable variations in the mid-IR spectrum, 
plus the metallicity of the star is unknown. The literature 
data on V42 (see Table appear to show a clear silicate fea- 
ture, which we have fitted using dusty. We have set the dust 
envelope's inner edge to be at 1200 K, near the tempera- 
ture limit for SiO conden sation under normal circumstances 
i|Nuth fc Ferguson! l2006| j. The fitted model has an optical 
depth of t v = 0.25, with the remaining input parameters 
staying the same as the model for V6. Assuming a metallic- 
ity of [Fe/H] = -1.62, this yields rates of M dust = 8.6 x 10" 11 
M yr" 1 and M tota i = 7.2 x 10" 7 M yr" 1 . 

The Gemini spectrum shows excess above the Rayleigh- 
Jeans tail that could be attributable to dust. The rise to- 
wards 8 fim could represent the SiO fundamental band in 
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Table 3. Mid-IR photometry of lo Cen V42. 



Phase 1 


Date 


Ref. 2 


Instru- 


A 3 


Flux 4 








ment 


(cm) 




0.00 


02 Jun 94 


a 


TIMMI 


10 


3.1 




Feb/Aug 97 


b 


ISOCAM 


9.6 


2.5 




Feb/ Aug 97 


b 


ISOCAM 


12 


1.8 


29.87 


02 Mar 06 


c 


MIPS 


24 


1.9 


30.03 


26 Mar 06 


c 


IRAC 


8 


1.6 


30.87 


29 Jul 07 


d 


IRS 


5-37 


2.3 


31.07 


27 Aug 07 


e 


T-ReCS 


8-13 


1.3 



Notes: (1) Phase in pulsation periods from /-band maximum, 
taken to be MJD 49506 after Dicke ns. Feast &: Lloyd Evans! 
l|l972l), with P = 148.64 day s (vL+ 00). (2) References: a - 
lOriglia. Ferraro. &: Fusi Peccil jl995h : b - O+02; c - B+08; d 
- G. Sloan, private communication; e - this work. (3) Ranges 
denote spectra. (4) Flux as a multiple of the photospheric con- 
tribution, extrapolated from spectral models; flux calculated 
using a F oc A - approximation for photometric observa- 
tions that are not at 10 (im. 

emission, but this is near the edge of the atmospheric win- 
dow and thus uncertain. The absence of a silicate feature 
in the spectrum is puzzling, given the clear feature appar- 
ently present in the literature broadband photometry (see 
Table [3] & Fig. I12|l . These historical estimates for the mid- 
IR brightness of V42 have been high. Previous estimates of 
its mass loss have suggested a mass-loss rate of 7-10 x 1CP 7 
Mq yr _1 , with a 200-350 K dust envelope, and it has been 
suggested that this is associated with a period of dust ejec- 
tion around 1925, possibly linked to a thermal pulse event 
(Origlia et al. 1995; O+02). This is similar, though cooler 
than, the mass-loss rate estimated from the dusty fit to our 
Spitzer and Gemini data. However, what is obvious from 
Table is the factor of 2.4 range in the 10-pm flux received 
from V42 since 1994. The obvious solution would be that 
this is connected with a periodic phenomenon linked to the 
pulsation cycle: there is not an obvious link, however, with 
pulsation cycle phase, but neither is there sufficient evidence 
to refute one. The amplitude of variation is also somewhat 
larger than the factor of 1.6 change in the L-band flux ob- 
served by Dickens et al. 1972. Note that all flux values in 
Table|3]show mid-IR excesses, strongly suggesting the pres- 
ence of circumstellar dust. 



Surprisingly for this M-type star, the emission from the 
Spitzer photometry and Gemini s pectrum are w ell-fit by an 
amorphous carbon wind (following Hanncr 1988h with a con- 
stant temperature of 600 K and a mass-loss rate of between 
1.7-4.3 x 10 -10 M yr" 1 in dust and 1.4-3.6 x 10~ 6 M© 
yr _1 in total. The variance here comes from the uncertainty 
in the absolute flux calibration of the spectrum; and the 
interpretation of carbon-rich dust also holds for the Spitzer 
IRS spectrum. Note that this figure is some 2-4 times greater 
than that derived from fitting silicate dust to the literature 
photometry. The carbonaceous nature of the wind is also 
puzzling, as this is not a bone fide carbon star (vL+07) . 

It is clear that our Gemini spectrum, and a similar 
Spitzer IRS spectrum taken four weeks prior (G. Sloan, pri- 
vate communication) show no silicate emission. What is in- 



teresting is that, although the spectrum does not change 
appreciably between the two epochs, the flux received from 
V42 decreases considerably, mirroring the wide range in the 
historic photometry measurements. We discuss possible rea- 
sons for this in Section \6. 5. 41 

The key consequence of the mid-IR variability observed 
in V42 is that it now becomes virtually impossible to ac- 
curately pin down a mass-loss rate to better than a factor 
of ~2 without simultaneous observations of the flux at ^ 4 
ium, where the dust component does not dominate, to define 
the stellar temperature and luminosity. 

With respect to the rest of our data, it would appear 
that the accuracy of the mass-loss rates of other stars may 
depend on their optical variability, plus the optical proper- 
ties of their circumstellar dust. This is a potential source of 
error in our analysis, but given that V6 and V42 are by far 
the most extreme pulsators, it would appear that V6 is the 
only other likely candidate for significant mid-IR variability. 
Our only comparison here is between the Spitzer photomet- 
ric data and the Gemini spectrum, but would appear to limit 
the effect to a level much less than the known errors present 
in our mass-loss analysis. 

5.3 Derivation of mass-loss rates for other stars 

5. 3. 1 Correction of systematic differences between Spitzer 
photometry and model atmospheres 

Fig- H3] shows the distributions of excess above the model at 
8 and 24 pm, respectively. Clearly, V6 and V42 are the most 
excessive stars, although there are a few others exhibiting 
excesses at both 8 and 24 /im. Those objects with flux ex- 
cesses at low luminosities are likely background galaxies or 
mis-identifications. 

There is an observable offset between the model at- 
mospheres provided by MARCS and the 8- and 24-^tm pho- 
tometry data for the warmer, less luminous stars (the giant 
branch is offset from unity towards the bottom of Fig. I13[) . 
This effect is only ~1% for the brightest stars, but increases 
to ~4% and ~10% for the fainter stars at 8 and 24 fim, re- 
spectively. The reason behind this is unclear, as from a mod- 
elling perspective one would assume that the Ray leigh- Jeans 
tail that has been used in extending the model atmospheres 
out to the limit of the 24-/im filter would be more accu- 
rate in a warmer atmosphere. A sili cate absorption fe ature 
caused by interstellar material (c.f. lEvans et al.ll2002l ) can 
also be ruled out, as it is only expected to affect our 8- and 
24-^tm photometry by about 0.5% in each case, based on a 
100 K interstellar medium at Av = 0.4 mag. A combination 
of SiO and H2O absorption in the region of the 24-^im filter 
bandpass may affect the cooler giants. It may be that this 
under-luminosity in the fainter stars is a result of imperfect 
data reduction of the Spitzer MIPS data near the sensitivity 
limit of the instrument. 

While this correction is not significant on the level of 
individual stars, it is important when considering the clus- 
ter as a whole: without a correction for this effect, the total 
mass-loss rate for the cluster would be grossly incorrect. In 
order to minimise the effects of statistical scatter in the pho- 
tometry, we have assumed that stars with less-than-expected 
8- and 24-fim fluxes have negative mass-loss rates. This is 
unphysical, but by making this assumption we can examine 
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Figure 13. Observed mid-IR excess (attributable to dust) at 8 
and 24 fim, derived from fitting to model atmospheres. Over- 
plotted is the line we have taken marking 'zero' dust mass loss 
from Eqs. (|6]l and Representative lcr error bars are shown to 
the right of the plots - stars which have mid-IR excess above the 
3ct level at 24 (im are marked by large squares. 



stars with low mass-loss rates in a statistical manner, with a 
minimum bias from random uncertainty in the photometry. 

To calculate the correction, we have assumed that the 
average star produces negligible dust, therefore will have no 
flux excess. We then split the data into luminosity bins of 0.1 
dex, and for each, calculate the average offset. These form 
a tight correlation, except at the lower and upper limits, 
where poor-quality data and significantly dusty atmospheres 
(respectively) deviate the median. Taking the intermediate 
points (representing 21 bins, containing 1624 stars at 8 /xm 
and 1207 stars at 24 ^m), we have computed a best-fit line 
and removed this from the data. These corrections are: 



-fobs 



Fob 



+ 0.01 log(L)- 0.035, 



8 /./■ Ill 



+ 0.025 log(L) -0.27. 



(6) 
(7) 



5.3.2 Dust Temperatures 

For the remainder of the stars we examine, we have only 
the 8- and 24-pm data from which we derive a mass-loss 



rate. At this point, we are forced to make some assump- 
tions about the mineralogy of the dust. We assume here for 
simpli city that all sta r s show 'astronomical' silicate dust, fol- 
lowing iDrahwT&Jjei (|l984l ). This thus means that our dust 
mass-loss rates will tend to under-estimate the total mass- 
loss rate of the star, as our wavelength coverage misses the 
12-^im AI2O3 and 13-^im MgA^CU bands s een at the start 
of the proposed dust condensation sequence l|Lebzelter et al.l 
2006). Carbon-rich dust is also not considered for the ma- 
jority of stars: if it is present, the assumption of silicate dust 
alone will again lead to an under-estimation of total mass- 
loss rate, as in the case of V42, by a factor of two to four. 

Our first priority was to derive an estimated dust tem- 
perature for the stars as, for a given optical depth, mass loss 
can vary significantly with dust temperature. In order to do 
this, we constructed a grid of dusty models at 4000 K, with 
similar dust properties as those used for V6, but with a pure 
'astronomical' silicate dust component and tv = 0.02 and a 
standard MRN distribution. The temperature of the inner 
edge of the dust envelope was varied from 200 to 1500 K, in 
steps of 100 K. 

The dependence of the calculated mass-loss rate on the 
assumed photospheric temperature of the star is relatively 
weak: for a constant IR excess and stellar luminosity, a com- 
parison of dusty models shows that M scales roughly as 
T 1 ' 25 in the region of 4000-4500 K, providing a reasonable 
(within 10%) fit to photospheric temperatures up to 6000 K. 
This factor takes into account the increase in wind velocity 
(from increased radiation pressure) that dusty calculates 
for higher-temperature stars. A velocity scaling relation of 
Voa oc J 11 / 3 can be shown to match the output from DUSTY 
to within 2% over temperatures from 3500 to 6000 K. We 
make these corrections in our final rates, but they do not 
add a considerable effect over the range of temperatures at 
which stars appear to be losing mass by dust-driven winds. 

Errors in the stellar temperature will have an impact 
here, though their primary influence is on the inferred lu- 
minosity of the star, since our dust mass-loss rates use pho- 
tometry taken from the Rayleigh- Jeans tail of the stellar 
spectra. This will add a 10% error in the mass-loss rates of 
the most-luminous non-variable stars, and up to 30% in the 
mass-loss rates of the most-luminous variable stars, due to 
their higher temperature uncertainties. Due to the stochas- 
tic nature of these errors, the effect on the sample as a whole 
will be greatly reduced from this. 

We can then derive the temperature of the dust, under 
the assumption that the underlying photospheric spectrum 
follows a Rayleigh- Jeans tail. The [8]-[24] colour of the ex- 
cess emission from the dust will then be directly indicative 
of the dust temperature, provided the chemical make-up of 
the dust does not change. For each model, we derive the 
[8] -[24] colour we would observe by convolving the dusty 
model spectrum with the Spitzer filter responses and tak- 
ing the ratio of the resulting fluxes. We then compare our 
observed [8]-[24] colour with the models and linearly inter- 
polate between them to find the dust temperature. 

Unfortunately, the excess flux from the dust compo- 
nent is typically of order of the errors within the photom- 
etry or less: only 30 stars - 27 cluster members and three 
non-members - show excess flux above the model SED at 
24 (jrn with a significance of over three standard deviations 
(for comparison, only four stars show a similar flux deficit, 
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all of which are under five standard deviations). The dust 
temperature is also very sensitive to the zero points of the 
photomet ry: the addi t ional factor of 1.015 in IRAC 4 (8 
[im) from iRieke et all (|2008h raises the temperature of the 
warmest dust surrounding V6 by 110 K and V42 by 340 K 
and even more dramatically alters the dust temperatures of 
all other stars. We were therefore only able to obtain only 
dust temperature estimates for these two stars, plus a hand- 
ful of others, listed in Table 0] LEID 40220 is omitted as, 
despite showing clear mid-IR excess, its photometric errors 
are still too large for a temperature to be determined. 

For many of these 30 stars, the [8]- [24] colour of the ex- 
cess emission is sufficiently high or low that it is outside the 
bounds of our test range (100-1500 K). We would not expect 
significant amounts of circumstellar dust hotter than 1500 
K, and we would expect the inner edge of the dust envelope 
to exceed 100 K unless the star had recently completed a pe- 
riod of mass loss and was now in (mass-losing) quiescence. 
The stars with very cool dust are all low-luminosity objects, 
suggesting they are blends with background galaxies (see 
Section [6.1.3p . Stars with warm dust may have emission in 
the IRAC 8-/im band. Several of the stars with evidence for 
warm dust are known variables. It is quite possible that, in 
the cases of V152, V138, V148 and V161, changes in the 
stellar luminosity may impact our derivation of mid-IR ex- 
cess and thus dust temperature. It therefore appears that 
most of the cluster's mass-losing stars have warm circum- 
stellar dust, with the most evolved, most IR-excessive stars 
(V6, V42 and possibly V17) showing probably slightly colder 
dust. 

It is noteworthy that we find much colder dust around 
V6 using only the Spitzer data points than we do when 
we also consider the Gemini T-ReCS spectrum. This is due 
primarily to the marginally different dust chemistry, which 
has a substantial effect on the derived dust temperature, but 
may be related to low-level variability of the star between the 
24- and 8-/j,m observations, or the presence of an additional, 
lower-temperature dust component. 

For the remainder of this analysis, we have assumed 
that the temperature of the i nner edge of the dust en- 
velope is the typical 1000 K l|Habind Il996t ) in all stars 
other than VI, V6, V17, V42, and LEIDs 39105 and 
45232, where we take the temperatures or lower limits as 
given in Table [4] In terms of mass-loss rates, this repre- 
sents a roughly average value out of the possible range, 
with a possible variation of ±80% between the coldest 
dust likely based on Table [4] and the hottest dust likely 

to form (;$ 1500 K, based on interferometric observa- 

i it -i | 1 

tions and c ondensati on temperat ures (Salpctcr 1977; Drainc 

198ll; iKozasa. Hasegawa fc Sekil 1 19841; iDanchi et all 1 1994 

Tuthill et al.ll2000l ; iGauba fc Parthasarathy||2004l )). 



5.3.3 Calculating the mass-loss rates 



the RGB/ AG B 



Using these assumed temperatures, we can derive the ex- 
pected excess flux at 8 /im by taking the ratio of the 8-/xm 
flux from each dusty model and dividing it by the flux from 
a similar dusty model with no dust. This is repeated for 
the 24-/im excess. Our mass-loss rate then scales similarly 
to Eqs. (gj and J5]), with: 



Table 4. Dust temperatures (based on pure 'astronomical' sil- 
icate dust) for stars with the highest computed dust mass-loss 
rates. Errors are on the basis of Spitzer photometry only and do 
not include the comparatively small uncertai nties in estimating 
stellar luminosity. Variable numbers are from lClementl lll997l) . 



LEID Variable 


Max. Dust 


Inner radius 


Notes 


number 


temperature 


of dust envelope 






(K) 


(R.) (AU) 





55114 
44277 
49123 
41455 
39105 
32029 
47153 
45232 
44262 
34041 
35250 
33062 
37098 
45169 
47135 
40015 
31159 
44457 
54135 
51019 



V152 
V138 
V148 

VI 
V161 

V42 
V2 
V17 
V6 



>1400 
>1400 
>1400 
>1400 
>801 
>685 
>595 
>556 

937 

715 

509 

320 



+207 
-143 
+ 109 
-88 
+ 177 
-138 
+48 
-48 

<706 
<371 
<328 
<240 
<186 
<100 
<100 
<100 



<23 
<50 
<26 
<27 
15 +» 
22 +» 
50 tti 
111 +« 



<8.8 
<8.8 
<8.4 
<10.7 

™ ±l:S 



23.2 
72.6 



+ 20.5 
-10.7 
+ 27.9 
-17.7 



Notes: (1) Post-AGB, proper motion membership probability 
3%, radial velocity 188 km s" 1 ; (2) proper motion and radial 
velocity non-member V825 Cen; (3) proper motion member- 
ship probability 0%, no radial velocity data, variable identified 
in vL+00 with P = 1.0253 days; (4) suspected SiO or other 
emission near 8 (im; (5) faint stars, suspected blends with 
background galaxies; (6) for silicate dust. 



Mi 



where: 



Mb 



200 



(JL) 

V200/ 



~ 1/2 ( 

Uo 4 J 



-Re, 



-Re = 



F°(8)+F°(24)- 
-F c d (8) +F c d (24) 



(8) 



(9) 



F° is the flux observed above the Rayleigh- Jeans tail in our 
SED fit at 8 and 24 fim, and F,f is a similar value, calculated 
for our dusty model at the relevant dust envelope temper- 
ature and scaled to the star's luminosity. The factor of 200 
in Eq. (JSJ) comes from the conversion from the dusty total 
mass-loss rate, which assumes a gas-to-dust ratio of 200, to 
a mass-loss rate purely for the dust component. 

Our computed mass-loss rates are listed in Table [SJ 
which excludes those stars identified as blends at 24 /xm 
in B+08. The correlation between mass-loss rate and lu- 
minosity is shown in Fig. 1141 which confirms that it is the 
most luminous objects that are losing observable quantities 
of dust, and that most of the dust mass loss is indeed com- 
ing from V6 and V42. We also note that V6 and V42 have 
very similar mass-loss rates in our analysis here, which only 
contains the Spitzer data, as in our analysis in Section T5. 3. 2 1 
which also includes our Gemini spectra and other data. This 
is despite the assumed dust temperatures being considerably 
different in both cases. 
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Table 5. List of mass-loss rates from corrected IR excesses, in 
order of decreasing total mass loss, excluding stars with known 
blending issues from B+08, stars with 24-fim excesses of > 3<r and 
stars dimmer than 500 Lq. Variable numbers are from Clement 
(1997). 
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Notes: (1) Values based solely on Spitzer photometry; (2) 
post-AGB star; (3) also shows significant (> 3<r) excess at 8 
fim; (4) shows 24-^tm excess below 3<r, but shows 8-fim excess 
above 3cr; (5) carbon star, see Section 16.1.21 (6) no metallicity 
known, we have used the cluster average, [Fe/H] = -1.62. 

Also noteworthy are the stars with apparent signifi- 
cant dust mass-loss rates at low luminosities (L ^$ 5OOL.0, 
M duBt ~ lO _11 M0yr _1 ) in Fig.Ql] Close inspection of these 
targets show that they are mostly warm stars (some on the 
HB) with very little or no 8-/xm excess, but considerable 24- 
/im excess. It appears most likely that these are cluster stars 
blended with background galaxies (see Section \6. 1. 3 1 . 



6 DISCUSSION 

6.1 Notable objects besides V6 and V42 

6.1.1 Post-AGB objects 

The cluster contains two confirmed post-AGB objects: 
Fehrenbach's Star (LEID 16018) and VI (LEID 32029). It 
is notable that, while Fehrenbach's Star does not appear 
to show a dusty circumstellar envelope, VI does. From the 
spectra of vL+07, Fehrenbach's Star also has very weak 
emission lines, but those of VI are strong. 

It is difficult to reconcile the similarities between the 
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Figure 14. Correlation between mass-loss rate and luminosity 
for confirmed proper-motion cluster members. Objects with 24- 
fira excess over 3<r are marked by crosses. 



dust shells in VI and the mass- losing AGB/RGB-tip stars 
(see Table given that VI has evolved off the AGB tip. 
It is worth noting that VI is almost four times more metal 
poor, so we may not necessarily expect the historic wind 
from VI to be the same as in V6. The proximity of VI 's 
dust shell to the star does raise the interesting possibility 
that the dust shell in VI is relatively stationary, supporting 
the low velocity rates given by dusty's scaling of v oc ip -1 / 2 . 

For Fehrenbach's Star, we can surmise that the absence 
of dust is due to either the dust having dispersed and become 
too cold to detect with our Spitzer data, or the dust having 
been destroyed. It does not appear to be present at any 
significant level in the lo nger- wavelength AKARI data from 
iMatsunaea et alJ <|2008h . 

Several other potential post-AGB stars have been high- 
lighted in our HRD (Fig. [3]). These are listed in Table H Of 
these, several show bad cross-identifications or poor quality 
2MASS data. As we expect such objects to scatter exces- 
sively in the HRD, this should not be surprising, and we an- 
ticipate these are the exceptions, rather than the standard 
in our dataset. In particular, LEID 34135 is identified as an 
RGB-a star, LEfD 40240 is a known blend in our 24-^tm 
data. The exact status of LEID 38171 remains unresolved: 
it shows a mid-IR deficit in our Spitzer 5.8- and 8-^im data, 
though it may be an early- AGB star with poor 2MASS pho- 
tometry. The remainder have SEDs visually consistent with 
their ascribed temperature and luminosities, although we 
note that, of these, only LEIDs 29064, 30020, 32015, 43105, 
46054 and 46265 have radial velocity measurements that 
confirm their cluster membership. These could potentially 
be post-early- AGB stars, which have left the AGB sequence 
before reaching the thermally-pulsating stages. None show 
a mid-IR excess at 24 /im above the photometric errors. 



6.1.2 Carbon stars 

Of the known carbon stars in the cluster (LEIDs 52030, 
41071, 32059, 14043, and 53019), only the most luminous 
of these, LEID 52030, shows significant dust mass loss, and 
even this is only a 4<r detection. The values we quote here are 
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Table 6. List of post-AGB candidates identified from the HRD 
(Fig-Oi showing luminosities, temperatures, proper motion mem- 
bership percentage probabilities (PM Mem.) from vL+00, radial 
velocities from vL+07 and individual notes. 
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> 18 953 
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Notes: (1) Poor or no 2MASS data; (2) blend in 2MASS; (3) 
bad cross-identification in 2MASS; (4) L ~ 190 Lq, probable 
early-AGB star with low mantle mass; (5) young post-early- 
AGB or very bright early-AGB star; (6) confirmed post-AGB 
star; (7) Fehrenbach's Star; (8) isolated object with very steep 
spectrum, nature unresolved; (9) likely early- RGB/ AGB ob- 
ject. 

for silicate dust. Carbonaceous dust reproduces the observed 
8-jLtm excess marginally better than silicate dust. However, 
any analysis of this star is uncertain as our photospheric 
models are for oxygen-rich stars. 

We estimate that the mass-loss rate for amorphous car- 
bon dust with the same grain parameters as our silicate 
dust model is roughly 50% higher, at ~ 4 x 10" 11 M Q yr _1 , 
though the error in the mass-loss rate itself is over 50%. 
It would appear at first sight that these stars are no more 
effective at producing dust than their M-type counterparts. 

6.1.3 Low-luminosity giant stars 

We find a number of stars below 500 Lq, both on the giant 
branches and the HB, which have significant IR excess. As 
they appear to occur with roughly equal prevalence across 
the HRD, we suggest that these are merely blends with back- 
ground galaxies (c.f. Section [5 . 3 . 2 [) . Fig. 1151 shows that the 
mid-IR-excessive stars nicely bifurcate into these two colour 
groups: we expect large 24-fim excesses but little 8-fj.m ex- 
cess ('cold dust') to most likely be blends with background 
galaxies (c.f. 2MASX J13272621-4746042 in B+08), while 
we expect fairly comparable 8- and 24-/xm excesses for stars 
with circumstellar dust. Note that V6 and V42 lie off the 
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Figure 15. Observed mid-IR excess (attributable to dust) at 8 
versus 24 fim, which corresponds (for given dust properties) to the 
inner temperature of the dust envelope. Potentially dusty stars 
and probable blends with background galaxies can be split by a 
line corresponding to a 600 K silicate dust envelope. Stars above 
500 Lq are shown as crosses, stars below 500 Lq as dots. 

diagram to the right, as does the low- luminosity RGB star 
LEID 31159, which lies just below the 600 K line. The low- 
luminosity stars with mid-IR excess also spatially follow the 
distribution of all stars, not just cluster members, and show 
no concentration towards the centre of the cluster, ruling 
out crowding in the cluster core as a primary factor for gen- 
erating these detections. 

In total, 37 stars with luminosities below 500 Lq show 
significant (> 3a) excess at 24 fim. Of these, however, 34 
have little or no excess at 8 /im. While we have made every 
effort to remove 24-/im blends from our sample, the glare 
from nearby stars up to 150 times brighter may prevent an 
accurate 24- ^im flux measurement being made in these cases. 
The remaining three stars are: LEID 31159, an 18 Lq RGB 
star; LEID 44457, a 49 Lq RGB star; and LEID 47135, an 
early-AGB star. We have no reason, a priori, to assume that 
any of these is likely to harbour circumstellar dust. 

6.1.4 Non-members 

Our SED modelling also includes a number of stars that have 
been confirmed as non-members, either by proper motion or 
radial velocity measurement. These include V2 (V825 Cen, 
LEID 34041), an emission-line M-giant with a fitted tem- 
perature of 3346 K and a strong IR excess . It wa s identified 
as a radial velocity non-member b ylFeastJ (|l965h and has a 
period of 236 days fjKaluznv et al.lll996T l. 

A fit to the IR excess in V2 is shown in Fig. 1161 with 
a 850 K blackbody. The temperature of this blackbody is 
comparable to the 715±gg 9 K we find in Section \5. 3. 21 This 
blackbody has a flux density scaling constant (L/R 2 ) that is 
ten times that of the underlying photosphere; though with- 
out a good distance measure, we cannot estimate its mass- 
loss rate. We cannot attain a good fit to the spectrum using 
oxygen-rich dust, due to the high flux in the IRAC 4.5- and 
5.8-/xm bands and the absence of any silicate-like feature 
(which would show up at 8 and 24 fj,m). Carbon-rich dust 
could provide a fit, though without a mid-IR spectrum of 
the star, we have no reason to believe it harbours such dust. 
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Figure 16. The SED for the radial velocity non-member V2 
(filled squares). Over-plotted is a 3500 K spectral model (solid 
line) and the same model with a 850 K blackbody (dashed line). 

The star LEID 37098 (VLR 0-123) was identified in 
vL+00 as a variable star with a period of around a day. It 
exhibits little excess at 8 fj,m, but considerable excess at 24 
fj,m. We estimate the mass-loss rate at 1.4x 10 -8 (d/kpc) -1,5 
Mq yr -1 , though the absence of 8-fim emission suggests that 
the star is most likely blended with a background galaxy. 

6.2 Total mass-loss rate of u Centauri 

Our calculation of the total mass-loss rate for u) Cen is sub- 
ject to broadly similar errors as discussed in Section T5. 2 1 

• 30% uncertainty in the dust mass-loss rate from internal 
errors in DUSTY; 

• 30% uncertainty in the dust mass-loss rate due to 
uncertain dust chemistry and grain size; 

• negligible uncertainty in the cluster's metallicity, which 
would affect the gas-to-dust ratio; 

• ^ 10% uncertainty due to errors in scaling the mass- 
loss rate of the stars to the photospheric temperature (see 
Section rO^j) ; 

• ^ 5% uncertainty due to error in the integrated flux 
density of the stars from fitting the SED, including errors in 
stellar temperature; 

• ~ 10% uncertainty due to error in the interstellar red- 
dening; 

• ~ 20% uncertainty due to error in the distance; 

• errors due to the corrections we have made to the pho- 
tometric excesses above. Though difficult to discern analyt- 
ically, we have altered the parameters of the constants in 
Eqs. © and {7} , and obtained a visual estimate of the error 
involved, finding it to be at most 30%; 

• uncertainty in the dust envelope temperature, although 
we do not expect the associated error in mass-loss rate to be 
as high as the 80% error in individual stars (Section 15. 3. 2[) , 
as this arises from (random) photon noise. However, since we 
do not have any sensible upper limits on dust temperature, 
bar the maximum dust condensation temperature around 
1500 K, the error involved may still be as high as ~30%. 

Added in quadrature, this gives a total error of 65%. 

For the remainder of this analysis, we have removed any 
objects thought to be 24-/xm blends on the basis of B+08. We 



include VI, which is listed as a proper-motion non-member 
in vL+00 (at a membership probability of 3%), but a ra- 
dial velocity member in vL+07. We do not include mass-loss 
rates from stars with L < 500 L©, as these appear not to 
be losing mass (Section l6.1.3|l . We also include the mass- loss 
rates derived from our Gemini spectra for V6 and V42, in 
preference to those found from Spitzer data alone. Integrat- 
ing the mass-loss rates of the remaining cluster members, 
we find the total stellar dust mass-loss rate within u> Cen, 
under the same assumptions listed in Section \5. 21 to be: 

Mdust = 1.3 ±g;§ xlO" 9 M Q yr" 1 . (10) 

Due to the incompleteness of our knowledge on precise dust 
composition, detailed in Section r5.3.2l we suggest that, in re- 
ality, this is a lower limit. It is notable that, of this value, V6 
and V42 therefore presently produce around half of the dust 
produced by the entire cluster, with V6 being the marginally 
larger contributor due to its higher metallicity. Also, the 
cluster's entire dust mass loss appears to come entirely from 
stars within a magnitude of the giant branch tip. This is in 
contrast to the case of 47 Tuc, where O+07 find that most 
of the dust mass loss occurs in the more numerous stars 
much further d own the giant bran ch. Even the shallow ex- 
posure of M15 l|Bover et al.ll2006l ) revealed dusty giants at 
least 1.5 magnitudes below the tip of the RGB. It is possible 
that O+07's result was affected by blending at longer wave- 
lengths in this visually more compact cluster, which would 
produce apparent excesses at these wavelengths. Observa- 
tions of other clusters are needed to confirm the extent of 
the dust mass loss. 

In our analysis, we have not considered gas mass loss by 
methods other than a dust-driven wind. The most important 
of these is chromospher ic mass loss. There is weak evidence 
l|Smith fc Dupreelll98Sl ; Mauas et al. (2006); MvL07) that 
the dust-driven wind begins following the extension of the 
stellar chromosphere by pulsation, although the exact se- 
quence of events is yet to be fully explained. For this reason, 
as well as those described in Section 15.3.21 we here derive 
only a lower limit to the total mass-loss rate from the clus- 
ter. Multiplying the dust mass-loss rate for each star by its 
gas-to-dust ratio (see Section f5 . 2 . 1 1) . we can again integrate 
the total mass-loss rate of all the dusty cluster stars, finding 
it to be: Mtotal Z 8 ±1 xlO -6 M yr" 1 . If we include the 
contribution fro m dustless (chromospherica lly-driven) winds 
as prescribed bv lSchroder fc Cuntzl (|2005l ). this lower limit 
rises to around Mtotal ~ 1-2±(!:| x 10~ 5 M Q yr -1 , depending 
on the point at which chromospheric mass loss is replaced 
by dust-driven and/or pulsation-induced mass loss. 

It is interesting to note here that the fraction of the total 
mass loss provided by V6 and V42 is somewhat less, due to 
the higher metallicity (hence presumed gas-to-dust ratio) of 
V6 and the inclusion of dustless winds. V42 now produces 
around 20% of the gas produced by dust-producing stars 
(not including mass loss from giants with chromospheres), 
but V6 is reduced to around 9% of this total. 

6.3 Mass loss along the giant branch 

6. 3. 1 The onset and evolution of dust formation 

While we have no clear measure in the data presented here 
of mass loss by means other than a dusty wind, we can 
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Figure 17. Computed total mass-loss rates and their dependence 
on log((/) (top panel), luminosity (centre panel) and stellar effec- 
tive temperature (lower panel). Large squares denote stars above 
L > 500 Lq with individually statistically-significant mass loss, 
based on a 3cr cut on 24- /im excess. 



place constraints on when dust production starts to become 
efficient and explore the parameters upon which it depends. 
Fig. [17] shows the variation of mass-loss rate with gravity, 
luminosity and stellar effective temperature. It is difficult 
to say what happens in denser, fainter, warmer stars due 
to lack of sensitivity in the photometry. Despite this, the 
fact that extreme dust mass loss occurs towards minimum 
gravity and maximum luminosity is obvious, although it is 



not clear which parameter is physically responsible. There 
are arguments for mass loss to be linked to any or all of 
these parameters. 

Surface gravity is a measure of escape velocity and of 
the ease with which dust grains can be accelerated away 
and lost from the star. It has previously been suggested 
that this is the primary determining parameter of mass loss 
i|Judge fc Stencellll99lf ). and we do indeed observe a smooth 
trend in Fig. Q7] 

Radiation pressure on the dust grains increases with 
increasing stellar luminosity (and likely temperature, de- 
pending on the grain opacities). Stars evolving along the 
giant branch increase monotonically in luminosity, but at 
the same time their photospheres cool as their radii increase, 
also in a monotonic fashion. As the stars in GCs have very 
similar masses, the differences in surface gravity amongst 
stars along the giant branch essentially depend on the ra- 
dius. Therefore, it is not surprising that if a correlation is 
found with any of the three parameters luminosity, temper- 
ature and surface gravity, that a similar correlation is found 
with the remaining parameters. We do observe an increas- 
ing trend in mass-loss rate with luminosity. The post-AGB 
star VI fits the luminosity relation better, possibly since its 
mass-loss rate was determined while it was still on the AGB 
at a similar luminosity, but very different surface gravity. 

Dust formation requires a sufficiently low temperature 
at still a sufficiently high pressure - a condition most easily 
met in the atmospheres of the coolest stars. At warmer tem- 
peratures, we see that no stars with T e g ^ 4500 K are pro- 
ducing dust significantly. Following the trend to cooler tem- 
peratures, however, we find that there is a clear transition to 
statistically-significant dust mass-loss rates for around half 
the stars below ~4300 K. By 3900 K, all stars seem to be pro- 
ducing dust at a substantial rate. The onset of dusty winds 
in lu Cen therefore seems to occur at a somewhat cooler and 
more luminous threshold than found in the more metal-poor 
M15 (Meszaros et al. 2008). 

The sudden transition at ~4400 K could correspond to 
a sudden transition to allow sig nificant dust producti on to 
occur in dense environments, as lWachter et all (2008) sug- 
gest occurs in carbon-rich stars. This may be a particularly 
important factor if the temperature gradient is small and 
pressure (i.e. density) gradient is high in th ese metal-poor 
stars, as suggested by iMattsson et all ([2008). The fact that 
there is a regime where some stars show significant mass- 
loss rates, while others show insignificant mass-loss rates 
could be due to episodic mass loss, where stars produce dust 
for increasing periods of time until they reach ~4000 K, at 
which point the continuous dust production may take place. 
The fraction of dust-producing stars at a given temperature 
would then yield the relative duration of the dust-producing 
phase. It may also be possible that the AGB stars (which 
typically have a lower surface gravity at a given surface tem- 
perature) can maintain a significant dusty wind already at 
higher surface temperatures than their RGB counterparts. 
The gas-to-dust ratio may also play a role: it is possible to 
have higher mass-loss rates at higher temperatures, but in 
winds containing fractionally less dust. 
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Figure 18. Comparison of our estimated total mass-loss rates 
to literature relations. Our data have been binned in units of 
0.05 log(LR/M) and are shown as boxes with heights of ±lcr. In- 
dividual stars with statistically-significant mass-loss rates have 
been shown in the top panel as filled squares. Compariso ns: 
iReimerj (l975h - top panel, dashed line; |Ori glia et al.l (|2007l l - 
top pan el, dotted line; ISchr oder &: Cunta feoOot) — middle panel, 
crosse s; lNieuwenhuiizen^T'de"*Jage?![l'990h - middle panel, plus 
signs; W achter et al.l ll2008h - bottom panel, crosses. 



6.3.2 Comparison with literature relations 

The real drivers behind mass loss are likely to be 
determined by a combination of these parameters. 
Such has been the supposition of previous empiri- 
cal determinations of mass-loss laws, such as those 



of iReimersl (Il975 


), |Nieuwenhuiizen & de Jagerl 


ISchroder & Cunta 1 ( 


20051). and O+07. 



Reimers' law is based on the following simple parame- 
terisation for chromospheric mass-loss: 
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where L*, R* and M* are the luminosity, radius and mass 
of the star in solar units, and 7/r, is a constant, generally 
assumed to be in the range 0.4 to 3 ( Stancliff e fe Jeffervl 
l2007h . but most commonly 0.5 ([Schroder fc Cunt j l20o'5n T 
which we assume here. The modified relation given by 
ISchroder fe Cunt j l|2005t ) is based on the energetics of mo- 
mentum transfer in a chromospheric outflow, namely: 
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where ys is suggested to b e (0.8 ± 0.1) x 10 
iNieuwenhuiizen fc de Jagerl l|l99lt give a parameterisation 
similar to that of Reimers, but fitted to the entire HRD: 
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(13) 

where ?7n = 9.63 x 10~ 15 M yr -1 . O+07 base their rela- 
tionship on a similar method to ours, using (K-[8\) colour 
excess to attribute dust mass- loss rates, and fit the following 
parameterisation: 
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with gravity (g*) and radius (R*) in solar units. As ijj and 
Vco vary among our stars, we here take C — 1. 

We compare the expected relations from these laws to 
our data in Fig. 1181 In this figure, we have assumed M = 0.8 
Mq . For the O+07 relation, we used ip = 2000, on the basis 
that u> Cen is an order of magnitude less metal-rich than 
47 Tuc, where O+07 assumed ifr = 200. The other relations 
have no de pendence on We als o compare the carbon-star 
relation of IWachter et all i|2008ft for the SMC (at [Fe/H] 
~ —0.7), namely: 
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Though Wachter et al. use C" = 2.34x 10~ 5 M yr -1 as their 
constant, based on their modelled scaling with metallicity, 
we would expect C 1 to be roughly 1 x 10 -5 M yr -1 for 
uj Cen, which we use in the following comparison of the 
aforementioned relations. 

We can scale our mass-loss rates to match the above re- 
lations given our assumption that M oc y/ip, with the solar- 
metallicity value being 200. This places V6 at ip = 2600 
and the majority of the cluster (probably including V42) at 
ip — 8300. We do not expect the gas-to-dust ratio to attain 
solar-metallicity values (i.e. ip ~ 200) in these stars, simply 
due to the lack o f condensable material which can form dust 
l|van Loonll200SF ). 

Fig. E] shows the relation from O+07 successfully re- 
produces the distribution of our dusty stars in the range 
7.5 x 10 5 ^ LR/M < 2.5 x 10 B times the solar value, sug- 
gesting that M oc y/ip does hold between 47 Tuc and u> 
Cen. However, the quality of the fit decreases dramatically 
once we include the dust-free stars and stars where we have a 
marginal detection of dust. The gas-to-dust ratios we require 
for V6 and V42 to be consistent with O+07 are if) m 800 and 



Giants in uj Centauri 21 



500, respectively. This may not be unreasonable, if V42 is 
part of the cluster's most metal-rich sub-population, though 
Origlia et al. themselves find several stars showing dusty 
emission above their own model, suggesting the relation does 
not describe these more unusual objects very accurately. We 
also note that our data do not reproduce the high mass-loss 
rates seen in O+07 below LR/M ss 7.5 x 10 5 times solar. 

The classic Reimers' relation is meant to describe mass 
loss in solar-metallicity environments where a chromosphere, 
not dust, is the primary wind driver. We should therefore not 
be especially surprised that it does not fit the upper giant 
branch of this metal-poor cluster well, where we expect the 
chromospheres to be disrupted and a dust-driven wind to be 
present. That it over-predicts the dusty mass-loss rate below 
LR/M ~ 7.5 x 10 times solar suggests that chromospheric 
mass loss may dominate here, while its under-estimate of 
mass-loss for greater values of LR/M confirms that a dusty 
wind is th e primary mass-losing me chanism here. The for- 
mula from ISchroder fc Cunta l|2005h is essentially a down- 
ward revision of Reimers' law, derived for similar environ- 
ments. This provides a similar conclusion, although the over- 
prediction around LR/M « 2 to 7.5 x 10 5 times solar is less 
severe. 

The iNieuwenhuiizen fc de Jagerl j 19901 ) relation pro- 
vides a slightly better match to the slope of the relation than 
the other laws, but grossly under-estimates the mass loss, re- 
quiring ip ~ 100 — 500 for the dustiest stars but ip < 100 for 
V6 and V42. It is, however, in agreement with the upper lim- 
its further down the giant branch. The values of ip required 
near the giant branch tip are unphysical, but perhaps this 
should not be surprising, given that Nieuwenhuijzen & de 
Jager themselves state that the relation does not work well 
for dusty stars. 

Empirical literature relations therefore do not seem to 
be able to reproduce the mass-loss rates we find. Without 
data on the gas-to-dust ratios in the wind, we cannot tell 
whether it is truly the mass-loss rate which changes, or 
merely the gas-to-dust ratio, which in turn influences our 
derived mass-loss rates. In Fig. 1191 we show two fits to the 
data, given by: 

/ T \ 8/3 

M = 3.3 x 10~ 22 J M yr-\ (17) 

and: 

/ T \ 5/4 

M = 4.5 x 10- 14 M yr" 1 . (18) 

These two fits show, respectively, the fit to all stars 
with LR/M > 80 000 and the fit to only those stars with 
significant mass-loss rates (defined, as above, by a > 3a 
flux excess at 24 fjxa). These two fits are markedly different, 
and highlight the necessity to treat all stars, not just those 
which are observed to lose mass, when creating relationships 
between mass loss and other parameters. 

We stress that Eq. (|17[1 is not a replacement for the the- 
oretical relations examined above, but simply a parametric 
fit to our data. Taking T oc R~ 2 , it is in re l ativel y good 
agreement with the model of IWachter et all (|2008T ). aside 
from a factor of ~10 in the constant, despite the fact that 
their relation is modelled on carbon-rich dust producers and 
not the oxygen-rich dust we have in the majority of uj Cen's 



1e-06 




30000 50000 100000 200000 400000 

LR/M (solar units) 



Figure 19. Correlation between total mass-loss rate and physical 
stellar parameters. The dashed line denotes the fit to the entire 
data above LR/M = 80 000, described in Eq. {TTJl; the dotted line 
fits only the stars with significant mass loss, described in Eq. 1181 . 
Individual stars with statistically-significant mass-loss rates have 
been shown as filled squares. 

stars. If we assume their almost-linear scaling of mass-loss 
rate with metallicity, our models have a product of ipVao 
that is about ten times too high. If we assume the scaling 
from dusty, the same is true for the product y/ipVoo. A de- 
crease in outflow velocity would be surprising, given how low 
the observed outflow velocities we calculate are (see Section 
ICj .3 . 3 p . A decrease to the value of ip by a factor of 10 would 
be more promising, though a factor of 100 would again bring 
ip to sub-solar values. 

6.3.3 Wind velocity 

Dusty prescribes by simple argument that the wind ve- 
locity (if the dust and gas are coupled) must scale as 
L 1 ^ 4 (ipp s )~ 1 ^ 2 , and also comes with the proviso that the 
mass-loss rates provided are only valid if the calculated 
wind velocity exceeds 5 km s" 1 . The main reason for this 
is that the turbulen ce in the wind is ex pected to be of or- 
der 1-2 km s _1 (e.g. ISchoier e t al. 2004), suggesting a wind 
that is highly fractal in nature, as large-scale turbulence 
caus es some parts of the wind to sink back towards the star 
(c.f. Ilcke. Frank fc Heskdll992l ). 

From our dusty models, we find the calculated veloc- 
ities for V6 and V42 are only 1.16 and 0.76 km s _1 , re- 
spectively (assuming [Fe/H]v42 = -1.62). A typical RGB- 
tip star, with L ~ 2000 Lq and [Fe/H] = -1.62, then has a 
wind velocity of ~1 km s _1 . Clearly this is much less than 
dusty's minimum valid wind velocity, and gives us some 
cause for concern. Observationally, very similar very-low- 
velocity w inds have bee n suggested in Galactic halo car- 
bon stars (|Mauronl l2008), as well as the likely more metal- 
rich Galactic M5-M6 Hie star L 2 Pup and the M8 IIIv 
star EP Aqr, which show hin ts of some asphericity or inho- 
mogeneity in their outflows l|Jura. Chen fc Plavchan| [2002: 
I Winters et"aT]|2007i ). 

Rotation and magnetic activity could also play a role 
in modifying mass loss. Perhaps the largest effect, however, 
is present in the stars near the AGB tip. Here, pulsation- 
induced shocks could increase the wind velocity consider- 
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ably beyond that which dusty predicts (c.f. lBowenlll988l ). 
Depending on the strength of the shocks, controlled by the 
piston amplitude of the pulsation, a wind speed of 10 km s _1 
could the oretically be obtained. This is similar to the situ- 
ations in iMattsson et all ()2008f ) . where the kinetic energy 
input (along with the amount of condensable carbon) con- 
trols the wind outflow: this is c orroborated by observational 
evidence dvan Loon et al . 2008). Shocks could facilitate dust 
formation, which would only occur when the pulsation be- 
comes strong enough to cause shocks capable of accelerating 
the wind. Given the highest mass-loss rates we find are in 
LPVs, this is perhaps not surprising. A side-effect of this 
is that, since mass- loss rate depends linearly on velocity, 
the mass-loss rates observed for these stars would be signif- 
icantly higher than the values we publish here. 



6.3.4 Evolutionary status of dusty stars 

A comparison of the relative abundance of stars on the RGB 
and AGB can provide evidence of the evolutionary stage of 
our dusty stars. To compute this ratio, we have used the rel- 
ative evolutionary rates (dt/dL, proportional to the number 
of stars in a given luminosity interval - dN/dL) along the 
Padova and BaSTI RGB and AGB isochroncs: 



Fagb(L) 
Ra/r(L) 



(dN/dL) j 



(dN/dL) agb + (d7V/dL) RGB ' 
(dN/dL) agb 



(dN/dL)t 



(19) 
(20) 



These ratios define, respectively, the fraction of stars on the 
AGB, and the ratio between AGB and RGB stars, at a given 
luminosity on the giant branch. In practice, these values 
are upper limits, signifying the fact that some stars may 
evolve off the AGB as a post-early- AGB or AGH-manque 
star before reaching the model's canonical AGB-tip (defined 
here by the last thermal pulse) . The Padova isochrones sug- 
gest Fagb declines from 30% to 20% between 100 and 2000 
Lg, while the BaSTI isochrones suggest the decline is from 
~19% to ~13%. We note that this percentage depends on 
the physics of convection used in the models, but that em- 
pirical constraints on convection mean that the physics used 
is unlikely to lead to significant alteration of our conclusions 
(Vandenberg et al. 2006). 

We have used the Padova isochrones as an upper limit 
for -Fagb in Fig. 1201 where we present the fraction of dusty 
stars as a function of statistical significance. Here, we can 
afford to relax our criteria for significance from our previ- 
ous value of 3(J, as we are dealing with a statistical sample, 
rather than individual stars, thus we expect errant data due 
to random noise will average out. 

This figure suggests that even from the start of mea- 
surable dust mass loss (~ 1300 Lq), at least some stars on 
the RGB must be producing dust, at a confidence level of 
~ 6<t. It also suggests that all stars which reach the RGB 
tip become dust producers. 

Both the Padova and BaSTI isochrones imply a reduc- 
tion in dN/dL by a factor of 2-7 as for luminosities above 
the RGB tip, suggesting we might see 9-30 stars in the 1000 
Lg above the RGB tip. In fact, we see none, with the pos- 
sible exception of the metal-rich V6. This includes all post- 
AGB objects, both confirmed and unconfirmed, which fur- 
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Figure 20. The fraction of dusty stars in the most evolved re- 
gions of the RGB+AGB, as a function of luminosity. The solid 
lines show the fraction of objects, above a given luminosity, with 
combined 8 + 24-^m excesses of at least one, two, three and six 
times the photometric errors, in order of increasing thickness. The 
dashed line (right axis) shows the total number of stars above that 
luminosity. The dotted line (left axis) demarcates the maximum 
likely fraction on the AGB: J F^QgdL. Dusty stars above this 
line therefore imply dust production on the RGB. 



ther suggests that the majority of AGB stars do not regain 
the luminosity of the RGB tip in u> Cen. 

We can see how this may occur if we consider the im- 
plied mass-loss efficiency on the RGB. Our isochrone fitting 
(particularly using the Dartmouth models) suggests that 
RGB/AGB stars initially have masses of 0.85 Mq, that typ- 
ically 0.20-0.25 Mq of that is lost on the RGB and around 
half that value on the AGB, leaving a remnant of around 0.5 
Mq . A significantly larger spread of RGB mass loss than this 
is required to explain the entire HB distribution. If mass loss 
on the AGB occurs at the same rate as on the RGB for a 
given luminosity, we can estimate that the AGB mass loss 
is given simply by the RGB mass loss, divided by the com- 
parative rate of evolution (and hence dN/dL) between the 
RGB and AGB, thus: 



Magb(L) 



Mrgb(L) 
Ra/r(L) 



(21) 



The value of R A /r(L) is 0.20-0.25 near the tip, where most 
mass loss takes place, implying a mass loss of 0.04-0.06 Mq 
on the AGB. In practise, the mass loss on the AGB is likely 
to be more than this, as it is easier to levitate mass off AGB 
stars as their mass (hence gravity) is less. This would suggest 
that many (perhaps most) stars either do not reach the AGB 
(AGH-manque stars) or 'peel off the AGB before they reach 
the tip. This further corroborates with the large number of 
low-luminosity post-early- AGB candidates we find in Sec- 
tion 16.1.11 though these may also be very-low-mantle-mass 
horizontal-branch stars on the transition to the AGB. This 
in turn again suggests that most dusty cluster members are 
actually .RGB-tip rather than AGB-tip stars. 



6.4 Comparisons between stellar groups 

We have taken the metallicity data available on to Ce n gi- 
ants from ICavrel de Strobel. Soubiran fc Ralitel (|200l[ ) and 
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Figure 21. HRD, using our luminosity and temperature val- 
ues, but sho wing interpolated literature spectros c opic m etal- 
licities from Cavrel de Strobel. Soubiran fc Ralitel d200lT) and 
I Johnson et al.l [120081 1. 
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Figure 22. HRD, showing stars with measured CN (top) and 
barium (bottom) line strengths from vL+07. Small points show 
stars with no measurements, circles show stars with measured line 
strengths, CN- and Ba-rich stars are shown as filled circles. 



Ijohnson et alj (|2008l ) and super-imposed them onto a HRD 
of our data. These are displayed in Fig. I21I In Fig. I22I we 
show a similar diagram, identifying those stars measured as 
CN- or Ba-rich in vL+07. 

As expected, the stars determined to be metal-rich oc- 
cupy a region much further to the cooler side of the giant 
branch. Of the stars with measured line strengths, the CN- 



rich stars do not appear to cluster to any one region more 
so than the CN-poor stars. This confirms that CN enrich- 
ment is primordial and not a result of the stars' evolution 
themselves. The Ba-rich stars are more likely to be found 
towards the cool side of the sequence, down to luminosi- 
ties below the HB, suggesting that these stars are probably 
metal-rich, and not associated with the AGB (c.f. vL+07). 
The strong temperature-dependence of the Ba A4554 line 
used to determine barium richness also suggests that the 
strong lines in the Ba-rich stars are due to the stars' low 
temperatures and higher metallicities, and that the relative 
abundances do not differ as much as is apparent in vL+07. 

Fig. [23] shows the variation of the same indices with 
dust mass-loss rate. Here we define two indices to measure 
the richness of CN and barium. These are taken relative to 
the demarcations for richness marked in Figs. 13 and 20(f) 
of vL+07, and are described thus: 



a = S3839 - 2.2(CH4300 -0.1), 
k = Ba4554 - 0.127(S3839 + 1.26), 



(22) 
(23) 



where S3839, CH4300 and Ba4554 are the line indices as de- 
scribed in vL+07. The values of a and k are positive for CN- 
and Ba-rich stars, respectively. Interestingly, there does not 
appear to be any significant variation of mass-loss rate with 
measured [Fe/H] values, unless one takes only those stars 
with highly significant excesses. However, there does appear 
to be a strong relation with barium-richness. This correla- 
tion is strengthened if we remove the special case of VI - the 
dusty post- AGB star. We also note that then, the four most- 
mass-losing stars are metal-rich, suggesting a weak correla- 
tion may still be present. The barium may simply trace the 
coolest and metal-richest subpopulations in uj Cen, which 
one would expect to be the most prolific dust producers. 

Using the above criteria, and setting a cut at 500 Lq, 
there is no obvious dependence of dust production on CN- 
richness: of the cluster members with observable mass loss 
(> 10~ 12 M yr" 1 ), four are CN-rich, and 53 CN-poor. This 
is comparable to the figures of 12 and 95 in the non-dusty 
sample. The only CN-rich star losing large amounts of mass 
is, perhaps unsurprisingly, the carbon star LEID 52030. The 
lack of CN-rich M-type giants has been noted in vL+07. 



6.5 Implications and further work 

6. 5. 1 Mass-loss rate and variability of metal-poor giants 
There are several implications that result from our analysis. 

• Mid-infrared excess, thus the inferred amount of dust, 
appears highly variable in pulsating stars on timescales of a 
month or less, much shorter than thought possible. Actual 
mass-loss rates may thus be considerably different from liter- 
ature values. The true mass-loss rate may not be attainable 
without mid-IR photometry and spectroscopy throughout 
the pulsation cycle. 

• An accurate measure of dust and gas mass loss requires 
knowing: how the gas-to-dust ratio varies both with evolu- 
tion and metallicity, whether gas is coupled to dust, and how 
low the outflow velocity can reasonably be to support strong 
mass loss. 

• Lack of correlation between our results and literature 
relations and studies (Sections l6.3.2l and l6.3.3[) highlights the 
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Figure 23. Dust mass-loss rates for stars with L > 500 Lq 
and their dependence on [Fe/H] (top panel), and the CN (cen- 
tre panel) and barium (lower panel) richness indicators a and k 
(see text). Large squares show stars with individually statistically- 
sigificant (> 3<r) mid-IR excesses. 



difficulty in obtaining empirical relations between mass loss 
and other stellar parameters in an inhomogeneous sample 
of stars. From our analysis of V42 (Section l5.2.2p and other 
stars (MvL07), it seems more likely that other non-LTE ef- 
fects, such as shocks and chromospheres, are significantly 
modifying the winds. 

• Our mid-IR spectrum and critical evaluation of other 
data of V42 suggest that the analysis of lHofner fe Andersen! 



(|2007l ) may be at least partially correct in assuming that 
some C-rich mass loss can occur in conditions that depart 
from LTE. This may provide t he opaci t y sou rce required to 
drive O-rich dust-driven winds IWoitkel (|2006l ) . 

6. 5. 2 The fate of u Cen's lost mass 

The fate of the lost mass is also interesting. The mass 
loss from the cluster must be stochastic at some level if 
it is being fed by only a very small number of stars, as 
these stars do not maintain a constant dust production 
throughout their evolution. Based on our observations, the 
current mass loss from the cluster should be around 1.2 
x 1CT 5 Mq yr _1 or more, depending on the contribution 
from chromospherically-driven winds, the gas-to-dust ratio 
in winds, and the gas-to-dust coupling efficiency. 

With this total mass-loss rate and assuming 48 Myr 
since the last Galactic plane crossing (when tidal friction 
is expected to strip away intra-cluster material) based on 
the proper motions from vL+00, we might expect <; 290 
Mq of gas and ~0.05 Mq of dust to have built up within 
the cluster. Our Spitzer observations (B+08) suggest the 
amoun t of dust in the cluster is < 1CP 4 Mq fc.f. lBover et al.l 
(2006), who detect 9 x 10~ 4 Mq of dust in M15), and data 
in the liter ature imply that < 2.8 Mq of neutral gas can 
be present l|Smith et al.lll990l ). The possibility remains that 
the intra-cluster medium (ICM) in ui Cen is in a hot, ionised 
state - the ICM in 47 Tuc contains a t least 5% ionised gas 
|Smith et al.lll990l ; iFreire et al] l200ll ). The fact that cold 
ICM in both gas and dust phase is observed in M15 suggests 
that neutral ICM is likely an existing phase in many GCs. 

On this basis, it would seem likely that neutral gas is 
being processed or removed from the cluster environment on 
timescales of J; 3 x 10 5 years, and du st within 1 x 10 s years. 
With a tidal radius of around 85 pc (|Harrislll996l ). it would 
take a wind 8 Myr to escape the cluster if it were escaping at 
10 km s~ , or 1.2 Myr if swept clean as a result of ram pres- 
sure by the hot Galactic Halo (u Cen has a s pace velocity 
of_69.1±5.6 km s _1 with respect to the Halo; lOkada et al.l 
120071 ). 

Another possibility is that the dust mass-loss rates (and 
thus also the gas mass-loss rates) we derive are too high due 
to their intrinsic variability, however we would expect this 
to average out through addition of errors in quadrature so 
that our value is only affected by this at the ~15% level, 
which mostly comes from the error in the mass loss of V42. 
Without knowing the long-term trend of this mass loss it is 
difficult to tell, but we can argue based on the sustainabil- 
ity of the mass-loss rate of V6, given the total mass loss we 
expect over its RGB or AGB evolution. Based on the total 
mass-loss rate we derive in Section [5.21 we find that, at the 
current rate, the star would lose the ~0.1 Mq required by 
evolutionary theory in 90 000 years. Given its position at 
the top of the giant branch, this does not appear too un- 
reasonable, so we can assume our global rates are relatively 
sound. 

We can also consider the dynamical effect of lost mass 
on the cluster. If a mass-loss rate of 1.2 x 10 -5 Mq yr _1 
is indeed representative, this eq uates to some ~0.46% of 
the c l uster's mass per gigayear (^M crritt. Mcvlar Tfc Mayor] 
Il997l : iD'Souza fc Rixl 120051 : | van de Ven et al] l2006h . If as- 
sumed constant in the past, the cluster will thus have lost 
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5-7% of its initial mass, with only small dynamical conse- 
quences. However, more massive stars lose a much larger 
fraction of their mass, and thus the cluster mass loss will 
have been more significant in the past - even when the clus- 
ter was already a Gyr old and intermediate-mass AGB stars 
were the main contributors to this mass loss. Cluster mass 
loss is therefore not a phenomenon limited exclusively to the 
first 10 Myr or so when exploding massive stars assured sig- 
nificant mass loss from the cluster's gravitational potential. 

6.5.3 Are globular clusters only producing dust 
episodically? 

Our analysis appears to show that the total mass loss from 
GCs is highly episodic. In u) Cen, V6 and V42 appear to 
be producing most of the dust mass loss in the cluster, and 
the total mass loss from the cluster appears dominated by 
a small number of stars. 

Mass loss in ui Cen would seem to be stochastically 
variable. The high number of low-luminosity post-(early)- 
AGB stars and candidates suggests many stars may leave 
the AGB early, and that there may be a considerable num- 
ber of AGB-manque stars, leaving a significantly depleted 
fraction of stars actually reaching the AGB tip. By impli- 
cation, the mass loss from other GCs of similar metallic- 
ity would have even more stochastically-variable mass-loss 
rates, as they are less populous and therefore have (numer- 
ically) even fewer stars exhibiting high mass-loss rates. 

Mass loss from a cluster may increase dramatically in 
the event that a single star undergoes the ejection of a plane- 
tary nebula. Several of these hav e now been o bserved in GCs: 
notably Pease 1 (K648) in M1 5 dPeasdl 19281), but also IRAS 



18333-2357 (GJJC1) in M22 (iGillett et al 



Pal6, and JaFu2 in NGC 6441 (|Jacobv et al 



1989), JaFul in 
1993 ). The fact 



that, out of all the clusters observed to date, only M15 shows 
convincin g evidence of hos ting an observable intra-cluster 
medium (|Bover et all 120061 ; Barmby et al. 2008), could be 
due to a random effect, rather than due to M15's partic- 
ularly low metallicity. This implies that the diffuse matter 
seen in M15 could merely be a diffusing planetary nebula 
that was ejected ^ 10 4 years ago, however its large mass 
0.3 Mq) would suggest that a single such event is unlikely to 
be the sole producer of this much material. Potentially this 
free gas and dust could also be the results of a stellar colli- 
sion, which by their nature are also episodic. Simulations in 
GCs fail to reproduce this leve l of mass loss in the case of 
collid ing main-sequence stars l|Lombardi. Rasio fc Shapirol 
1996), however a collision between a red giant and a neu- 
tron star may have the potential to eject the entire ~ 0.3 
M© stellar mantle (|Rasio fc Shapirol [l99lh . Such collisions 
have been suggested as a method o f clearing GCs of ISM 
jllmbreit. Chatteriee fc Rasiol I2008T I . The true fate of the 
lost mass, and the precise origin of the cloud in M15, re- 
mains at this time unknown. 



6.5.4 The nature of V42's mass loss 



It is clear from Section 15.2.21 that V42, with its mass- 
loss rate of 1.4-3.6 xlO -6 Mq yr _1 , stands out from the 
rest of the group. One possibility for this is that this rep- 
resents a 'superwind' phase akin t o that seen in more- 
massive long-period variables (e.g. Ivan Loon et al.l 1 19991 ; 



ISchroder. Winters fc Sedlmavrl 19991 ). Here, the mass loss in- 
creases dramatically shortly before the star makes the tran- 
sition from AGB-tip star to post-AGB star. If this were the 
case, we would expect V42 to denote the AGB-tip in Fig. [3j 
As we are unable to separate the RGB from the AGB, we are 
unable to tell whether this is the the RGB tip may 

be more luminous. V42 is also clearly less luminous than 
V6 both in the data we have used to determine its luminos- 
ity, and the variability studies of Dickens et al. (1972). We 
should also note that the luminosity of V42 is likely vari- 
able on timescales significantly longer than the pulsation 
timescale - it is possible that V42 has recently experienced 
a thermal pulse, as suggested by Origlia et al. (1995), and 
that it is currently in an associated luminosity minimum. 

The position of V42 on the HRD (Fig. [3jl also suggests a 
second possibility - that the star is part of one of the metal- 
rich sub-populations, like V6. If it is 0.5 dex more metal- 
rich, this would decrease the calculated outflow velocity of 
its wind and hence its dust mass-loss rate by 0.25 dex and 
decrease its gas mass-loss rate by 0.75 dex to better match 
V6 in Fig. \M 

The possible presence of carbon-rich dust highlighted 
in Section 15.2.21 is surprising, given the oxygen-rich nature 
of the star. We put forward three suggestions to explain 
this behaviour: firstly, that a large variation in dust pro- 
duction is producing significant variations in the amount 
of dust in the wind; secondly, that luminosity changes in 
the star are altering the amount of light reprocessed by the 
dusty wind; or finally, that the grain size varies sufficiently 
that the silicate feature disappears periodically. In the first 
case, it can be argued that the criteria for dust formation to 
occur are only met during part of the pulsation cycle, lead- 
ing t o episodic dust formation under non-LTE conditions 
(e.g. iHofner fc Andersen! I2007T I . This may lead to carbon- 
rich dust at some epochs, with oxygen-rich dust forming 
at others, which we propose as a possible solution to the 
disparity between the mid-IR spectra (suggesting carbon- 
rich dust) and accumulated literature photometry (suggest- 
ing oxygen-rich dust). 

Under the second scenario, the rapid, large-scale 
changes seen in the mid-IR amplitude of V42 are due to 
the variation of light from the pulsation cycle: during peri- 
ods when the star is warm, there is more optical emission 
to be re-processed by the dusty wind, which is much more 
opaque at these wavelengths. During periods when the star 
is cool, the optical emission is suppressed, the re-processing 
of optical radiation is decreased, and the mid-IR emission 
drops accordingly. 

It would, at face value, appear unlikely that the grain 
size in V42 could become sufficient (i.e. a few tens of mi- 
crons) to completely mask the silicate feature, due to the 
low metallicity and thus the infrequency of particle collisions 
that would lead to the aggregation of large grains. However, 
we note that it is not presently possible to observationally 
distinguish between these three scenarios. 

Episodic and/or variable mass loss near the RGB /AGB 
tip i s not a new concept: previous studies, in cluding O+02 
and iMeszaros. Dupree fc Szentgyorgyil (|2008l ) , have argued 
on observationally-based grounds that episodic mass loss is 
occuring. Applying this principle to other stars, the conse- 
quences of this could be quite far ranging: historical mass- 
loss rates have been assumed to be accurate probes of the 
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mass lost over many years, yet many observed stars have 
very strong optical variations. It is thus important that, 
when calculating dusty mass-loss rates from such stars, that 
investigators have a good estimation of the stellar luminosity 
and temperature at the time of measurement. 



7 CONCLUSIONS 

We have here presented stellar parameters derived from 
spectral energy distribution fitting to stars over two orders 
of magnitudes in luminosity down the RGB and AGB of lu 
Centauri, creating a physical HRD of luminosity versus tem- 
perature, outlining the RGB, HB and AGB in great detail 
and accuracy, and identifying several new post-AGB star 
candidates. From isochrone and spectral energy distribution 
fitting, we estimate the following parameters for the cluster: 

• Distance: d = 4850±200 (statistical) ±120 (systematic) 
pc; 

• Reddening: E(B - V) = 0.08 ± 0.02 mag (statistical) 
±0.02 mag (systematic) ±0.02 mag (differential); 

• Total dust mass loss: M dust = 9 ±1 xl0~ 10 M Q yr _1 ; 

• Total gas mass loss: M gas ~ 1.2 x 10" 5 M yr _1 (~ 2/3 
from dusty winds, ~ 1/3 from chromospheric mass-loss); 

• Timescale to clear the intra-cluster medium from clus- 
ter: < 10 5 yr; 

under the following assumptions: 

• gas-to-dust ratio, ip oc 10~ [Fc/H ' ; 

• wind velocity, v oc L 1 ^ 4 ^ 1 ^ 2 ; 

• absence of substantial mid-IR variability; 

• constant dust chemistry along the giant branch; 

• inner dust envelope temperature of 1000 K; 

• our low wind velocities (~ 1 km s _1 ) are accurate and 
still produce an accurate mass-loss rate. 

We derive inner dust envelope temperatures for a handful of 
stars (including V6 and V42), suggesting they are typically 
> 600 K. 

We show that V6 and V42 contribute ~25% of the clus- 
ter's dust production. V42 may be variable at mid-infared 
wavelengths on timescales of as little as a few weeks, sug- 
gesting equally variable dust production, meaning a single 
estimate of mass-loss rate and wind conditions may not ac- 
curately reflect the long-term conditions in this and simi- 
lar stars. V42, an oxygen-rich star, may also be producing 
carbon-rich dust. 

Dusty mass loss appears to start suddenly at a threshold 
of ~4400 K, or ~ 1000 Lg, for a significant fraction of stars. 
There appear to be additional star-to-star differences, sug- 
gesting that other factors also influence the mass-loss rate. 
Empirical literature relations do not accurately reflect the 
dust production rates we measure in lu Cen. 

By comparing our stellar parameters with optical line 
strengths, we have deduced that the Ba-rich stars in lu Cen 
belong to the more metal-rich populations and not to the 
AGB. These stars are characterised by dustier winds. 

Mass loss along the RGB appears quite efficient with 
typically 0.20-0.25 M lost on the RGB and >0.05 M on 
the AGB, possibly leading to large numbers of AGB- manque 
and post-early-AGB stars. Dust production in lu Cen ap- 



pears not confined to the AGB, with most dust producing 
stars in the cluster near the RGB-tip. 
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